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PREFACE 
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1.0  INTRODUCTION 
1.1  BACKGROUND  AND  OBJECTIVE 


A problem  which  has  long  been  of  importance  to  the  experimentalist 
is  the  determination  of  radially  dependent  physical  properties  in 
axisymmetric  sources  from  measurements  of  integrated  properties  along 
a line  of  sight.  One  important  example  of  this  type  of  problem  is 
the  calculation  of  the  radial  distribution  of  emission  coefficients 
from  observations  of  the  radiance  (radiated  power  per  unit  area  per 
unit  solid  angle)  from  a cylindrically  symmetric  radiating  source.  When 
the  source  is  optically  thin,  the  solution  to  a problem  of  this  type  is 
usually  reached  by  the  use  of  the  Abel  transform  (Ref.  1).  The  appli- 
cation of  the  Abel  transform  in  this  instance  requires  the  emission 
coefficients  to  have  cylindrical  symmetry.  With  the  cylindrical  symmetry 
assumption,  the  measured  radiance,  I(x),  can  be  written  (Fig.  1)  as 


I (x)  = 2 J 

X 


^ f (r)  r dr 

. 2 2 'A 

(r  — x ) 


(D 


Here  I(x)  is  the  measured  radiance  (radiated  power  per  unit  area  per  unit 
solid  angle) , which  is  a function  of  the  displacement  x,  r is  the  local 
radius,  R is  the  maximum  radius  of  the  source,  and  e(r)  is  the  emission 
coefficient  (radiated  power  per  unit  volume  per  unit  solid  angle).  Note 
that  I at  each  x is  the  result  of  integrating  the  emission  coefficient 
e across  the  extent  of  the  source.  The  factor  2 arises  because  of  the 
cylindrical  symmetry  assumption.  Equation  (1)  is  a form  of  the  Abel 
integral  equation  (Ref.  2). 


The  use  of  Abel's  transformation  yields 


( (r) 


_1_  f " dl/dx 

11  . 2 2 !4 

Kx  — r ) 


d x 


(2) 
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Equation  (2)  is  the  inversion  equation  which  gives  the  radial  emission 
coefficient  in  terms  of  the  geometry  of  the  source  and  the  measured 
radiance  distribution.  Details  of  the  development  of  Eqs.  (1)  and  (2) 
are  found,  for  example,  in  Ref.  1. 


Figure  1.  Geometry  of  the  axisymmetric  source. 

A number  of  methods  have  been  applied  to  the  solution  of  Eqs.  (1) 
and  (2)  to  obtain  the  radial  distribution  of  the  unknown  emission  coeffi- 
cient e(r).  An  early  approach  documented  by  Pearce  (Ref.  3)  depends 
upon  finding  the  areas  of  homogeneous  zones  in  the  plasma  and  replacing 

the  integral  by  a summation  over  the  zones.  Nestor  and  Olsen  (Ref.  4) 

2 2 

make  a transformation  of  variables  such  that  r = v and  x = u and  assume 
that  I' (u)  is  constant  over  each  small  interval.  The  subsequent  series 
of  integrals  may  then  be  evaluated,  and  Eq.  (2)  is  approximated  by  a 
sum.  Bockasten  (Ref.  5)  fitted  third-degree  polynomials  exactly  to  the 
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data  points  and  approximated  Eq . (2)  with  a sum.  Numerical  methods  such 
as  those  described  in  Refs.  3 through  5,  which  do  not  use  smoothing  of 
the  raw  data,  have  an  intrinsic  disadvantage  in  that  small  errors  in  the 
radiance  can  lead  to  fairly  large  errors  in  the  emission  coefficient 
because  of  inaccuracies  in  the  derivative  of  I(x). 

Some  investigators  have  used  smoothing  techniques  utilizing  curve 
fitting  or  other  mathematical  approximations  in  an  attempt  to  reduce  the 
effect  of  experimental  error  in  the  data.  Freeman  and  Katz  (Ref.  6) 
least  squares  fitted  a single  curve  to  the  raw  data.  Barr  (Ref.  7)  used 
least-squares  polynomials  to  determine  the  best  fit  of  the  data  over  a 
number  of  intervals  centered  about  each  data  point.  Birkebak  and  Cremers 
(Ref.  8)  used  a method  similar  to  that  of  Barr.  The  data  were  least-squares 
curve  fitted  to  polynomials  over  a number  of  data  intervals.  Dooley  and 
McGregor  (Ref.  9)  directly  applied  the  integral  in  Eq.  (2)  by  using  the 
experimental  data  and  a coordinate  transformation  of  the  integrand  to 
provide  a means  of  numerically  solving  the  integral. 

Comparisons  of  various  techniques  (Ref.  8)  indicated  that,  in  general, 
smoothing  techniques  yield  better  final  results  than  other  methods,  parti- 
cularly when  there  is  appreciable  scatter  in  the  data. 

A major  problem  with  any  method  of  solving  the  problem  described  by 
Eq.  (2)  is  that  of  determining  the  effect  of  experimental  error  in  the 
measured  radiance  values  upon  the  resultant  values  of  the  emission 
coefficient.  In  an  effort  to  obtain  error  propagation  information,  some 
investigators  have  applied  their  methods  of  solution  to  simple  problems 
in  which  analytic  solutions  could  be  determined  (Refs.  8 and  10).  By 
using  scattered  and  unscattered  input  data  (i.e.,  radiance  values)  and 
by  comparing  the  resultant  emission  coefficient  values  with  their 
analytic  values,  the  investigators  were  able  to  obtain  empirical  infor- 
mation on  the  error  propagation  characteristics  of  the  techniques  as 
applied  to  particular  test  problems.  However,  at  present,  no  general 
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analytic  method  has  been  developed  which,  after  a prior  analysis  of  the 
input  data,  follows  the  propagation  of  experimental  error  throughout  the 
numerical  steps  of  the  specific  technique.  It  is  the  purpose  of  this 
investigation  to  describe  a method  for  smoothing  the  data  for  inversion 
and  to  develop  a concomitant  error  propagation  analysis. 

1.2  CRITERIA  FOR  CURVE  FIT  CHOICES 

Before  proceeding  to  the  description  of  the  specific  technique,  it 
would  be  useful  to  consider  some  of  fhe  more  subtle  consequences  of  curve 
fitting  in  order  to  define  a criterion  by  which  the  curve  fit  choices  used 
in  the  work  reported  herein  were  made. 

Experimental  data  are  subject  to  random  uncertainties  and  are 
generated  by  a physical  phenomenon  for  which  one  does  not  necessarily 
know  the  functional  form.  Indeed,  specification  of  the  functional  form, 
or  an  "adequate"  approximation  thereof,  is  generally  the  objective  of 
the  analysis.  Because  of  the  experimental  scatter  in  the  data,  one  may 
be  uncertain  of  its  correct  functional  form.  The  usual  approach  in  such 
situations  is  to  curve  fit  the  data  to  some  smooth  functional  form  either 
by  analytical  or  graphical  techniques.  When  the  analytic  approach  is  used, 
the  data  are  often  fit  to  a specific  polynomial  in  the  independent  variable. 
However,  as  is  often  readily  apparent,  the  first  choice  of  polynomial  often 
cannot  do  an  "adequate"  job  of  fitting  the  data,  and  other  polynomials  are 
tried.  This  failure  of  fit  may  be  caused  by  the  fact  that  the  data  are 
simply  not  expressible  by  a polynomial  function.  Similarly,  for  the  same 
reason,  there  are  often  several  choices  for  polynomials  which  appear  to 
"adequately"  fit  the  data.  Yet,  each  polynomial  is  different  and  will 
yield  somewhat  different  results.  Hence  the  results  of  a smoothing  process 
cannot  be  considered  unique  insofar  as  the  underlying  physical  phenomena 
are  concerned,  and  there  can  be  a large  number  of  possible  solutions. 

(Note  that,  for  a given  set  of  data  and  a specific  polynomial  form,  the 
numerical  solution  is  unique;  it  is  the  resultant,  derived  description  of 
the  underlying  physical  phenomena  which  is  not  unique). 
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The  observation  that  it  is  generally  unreasonable  to  expect  a 
polynomial  to  provide  an  adequate  representation  of  some  physical  phenome- 
non over  a wide  range  of  the  independent  variable  prompted  the  approach 
described  herein.  The  data  are  divided  into  several  intervals,  and  a 
different  polynomial  is  assumed  valid  over  each  interval.  Since  physical 
phenomena  are  generally  smooth,  the  polynomials  are  constrained  to  be 
smooth  and  continuous  at  interval  boundaries.  Since  the  data  are  expected 
to  be  randomly  scattered,  a final  least  squares  constraint  is  imposed  to 
provide  smoothing  capabilities. 

If  the  experimental  error  is  propagated  through  a curve  fitting  process, 
one  can  obtain  an  estimate  for  the  error,  or  uncertainty  in  the  results. 

Even  though  the  correct  functional  form  may  not  be  that  chosen  for  fitting, 
the  errors  are'  propagated  as  though  the  chosen  function  is  the  correct 
one.  If  the  chosen  function  provides  an  "adequate"  representation  of  the 
raw  data,  it  is  expected  that  the  derived  results  provide  an  "adequate" 
representation  of  the  true  physical  phenomenon.  Thus  the  propagated  errors 
induced  by  the  data  uncertainties  are  an  "adequate"  representation  of  the 
uncertainties  in  interpretation  of  the  physical  phenomenon.  For  the 
application  described  herein,  when  several  functions  appear  to  give 
"adequate"  curve  fits,  their  results  and  propagated  uncertainties  agreed 
well  with  each  other.  This  observation  makes  a strong  favorable  argument 
that  the  results  are  in  fact  descriptive  of  the  underlying  physical  phe- 
nomenon. However,  it  must  be  remembered  that  the  propagated  uncertainties 
describe  the  uncertainties  in  the  data  fitting  the  chosen  form  and  not  the 
uncertainties  in  the  chosen  function  fitting  the  physical  phenomenon. 

Since,  for  a given  set  of  data,  there  may  be  a wide  choice  of  possible 
curve  fits,  the  choice  of  which  curve  fit  to  use  is  highly  subjective  and 
generally  must  be  left  to  the  investigator.  The  criteria  for  determining 
which  result  to  use  are  as  follows: 
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1.  Choose  the  set  of  functions  which  perforins  an  acceptable 
fit  to  the  data,  based  on  the  observed  error  propagation 
characteristics. 

2.  From  the  set  of  functions  chosen  by  (1),  pick  the 
function  which  provides  the  best  error  propagation 
characteristics  for  the  problem. 

Although  the  foregoing  discussion  may  appear  obvious  to  the  experienced 
data  analyst,  the  underlying  theme  has  been  that  there  can  be  no  procedure 
or  set  of  conditions  which  will  unequivocably  define  the  correct  curve 
fits  in  all  cases.  The  error  propagation  analysis  provides  only  additional 
objective  data  for  what  is  finally  a subjective  decision. 

1.3  TECHNICAL  APPROACH 

The  numerical  approach  is  based  on  a curve-fitting  method  which  uti- 
lizes a least-squares  polynomial  spline  fit  technique  (Ref.  11)  to 
smooth  the  data.  The  Abel  inversion  is  effected  by  dividing  the  raw 
data  into  several  small  intervals  and  obtaining  a least-squares  polyno- 
mial curve  fit  to  the  data  in  each  interval.  Since  the  integral  is  a 
linear  operator,  the  emission  coefficient  integral,  Eq.  (2),  is  expressed 
as  the  sum  of  several  integrals;  each  integral  applies  over  a different 
interval  of  data,  and  in  each  case  the  radiance  data  are  expressed  by  a 
different  polynomial.  The  emission  coefficient  is  expressed  as 


1 

#*Z  J-l  dIj_l  / dx 

n 

__L  V 

r7>  i|i/d"  j. 

(3) 

77 

1 U A 

J 2 2 M 

r lx  — r ) 

77  2~l 

i=i 

J ,2  2 a 

7 • i ~ r ) 

where  the  subscript  i identifies  the  particular  interval  for  evaluation, 
and  denote  the  endpoints  of  the  interval,  n is  the  number  of 
intervals  required  to  span  the  data,  and  in  all  cases  r is  less  than 
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For  the  values  of  the  emission  coefficients  to  be  physically  correct, 
the  polynomials  must  be  constrained  to  be  smooth  and  single  valued  at 
interval  boundaries  and  also  must  satisfy  the  assumption  of  cylindrical 
symmetry.  The  cylindrical  symmetry  assumption  is  easily  accounted  for 
by  using  an  even  function  of  the  form 


2 

I|(x)  = u j | + a]7  x“ 


4 6 

a,3x  + ai4x 


(4) 


where  the  subscript  i denotes  the  ith  interval.  The  functions  in  adja- 
cent intervals  are  designed  so  that  their  ordinate,  slope,  and  second 
derivatives  have  the  same  value  at  the  endpoints  of  adjacent  intervals. 
Thus,  the  coefficients  of  a sixth-degree,  even  polynomial  over  each  inter- 
val of  data  are  determined  such  that  a best  fit  in  the  least-squares 
sense  over  the  entire  set  of  data  is  obtained  with  the  condition  that 
the  polynomials  and  their  first  derivatives  are  smooth  and  their  second 
derivates  are  continuous  at  the  interval  boundaries.  With  I(x)  repre- 
sented by  a series  of  polynomials  in  the  form  of  Eq.  (4),  the  integrals 
become  expressible  by  direct  integration  in  closed  form  over  each  inter- 
val of  the  spline  fit. 


Another  constraint  upon  the  functional  representation  of  the 
physical  data  which  is  imposed  by  the  development  of  Eq.  (2)  is  that  it 
must  have  zero  slope  and  ordinate  at  the  boundary  of  the  cylinder.  This 
constraint  is  imposed  artifically  by  obtaining  a new  curve  over  the  last 
interval  after  the  other  intervals  have  been  fit  with  the  proper  constraints. 
The  constraints  of  continuous  first  derivative  and  smooth  ordinate  are 
maintained  at  the  interval  boundary  for  this  artifical  curve.  The  con- 
struction of  the  artifical  curve  reduces  the  accuracy  of  the  curve  fit 
intensity  values  in  the  last  interval.  However,  the  curve  fit  values  of 
all  the  other  intervals  are  unaffected,  and  in  general  the  data  in  the  last 
interval  are  very  near  to  zero,  quite  inaccurate,  and  insignificant  with 
respect  to  the  rest  of  the  data.  When  the  results  of  the  inversion  in  the 
last  interval  are  important,  the  method  described  herein  must  be  applied 
cautiously.  A computer  program  to  perform  the  analysis  developed  herein 
is  described  in  Appendix  A. 
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2.0  NUMERICAL  INVERSION  TECHNIQUE  AND 
ERROR  ANALYSIS 


2.1  NUMERICAL  TECHNIQUE 


As  was  indicated  in  Section  1.0,  the  problem  which  prompted  the  present 
investigation  was  the  development  of  a method  to  analyze  the  effect  of 
experimental  error  in  observed  radiance  data  upon  the  emission  coeffi- 
cient values  obtained  by  the  solution  to  the  Abel  inversion  problem.  A 
measure  of  the  experimental  error  in  the  radiance  data  is  provided  by  the 
standard  deviations  of  the  radiance  values  which  are  determined  from  the 
experimental  data. 

From  statistics  theory,  it  is  known  that  if  a matrix  vector 
equation,  A = MB,  can  be  written  where  M is  a transformation  matrix 
relating  the  vectors  A and  B,  then  a relationship  between  the  standard 
deviations  of  the  elements  of  vectors  A and  B can  be  derived.  The  first 
step  in  the  error  analysis  is  the  development  of  the  matrix-vector  equa- 
tion, E = FY,  where  Y is  a vector  containing  the  radiance  data  values,  and 
F is  a matrix  relating  the  vectors  Y and  E. 


Since  the  data  are  to  be  divided  into  intervals,  denote  the  displace- 
ment values  which  are  endpoints  for  the  interval  by  ^ and  and 
k = 1,...,  n.  Note  that  for  n intervals  there  are  (n  + 1 ) z values, 
j = 0,...,  n.  The  polynomial  in  each  interval  can  be  written  as 


4 

Pk  U)  = aki  x2  ; k = 1,  2,  . . . , n 

i=  1 


(5) 


zk_i  < x < Zk 


where  n is  the  number  of  intervals  into  which  the  data  points  have  been 
separated.  The  coefficients  a^  will  in  general  be  different  for  the 
polynomials  representing  each  interval. 
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Since  the  raw  data  in  each  interval  are  to  be  least-squares  curve 
fit  to  polynomials  of  the  form  of  Eq.  (5),  an  expression  for  the  sum 
of  the  squares  of  the  deviations  of  the  curve  fit  values  from  the  actual 
(input)  radiance  values  is  needed  for  each  interval.  For  this  purpose,  let 


nik 

Sk  ( ak  ]•  ak2-  ak3'  ak4)  =2^  [i(xj)  - Pk  (\j)]  :k=l,2 n 

i = 1 

Z k-  1 ^ x:  ^ ^k 


(6) 


where  n is  the  number  of  intervals,  m,  is  the  number  of  data  points  in  the 

th  ^ 

k interval,  and  the  displacements  are  numbered  independently  in 

each  interval.  The  value  of  in  each  interval  must  now  be  minimized 

subject  to  the  constraints: 


^k,l  (akl 

’ uk2'  ak3’  ak4 ^ 

= w 

- pk+ 

l«k> 

= 0 

?k.2  (ak]’ 

uk2’  ak3’  ak4^ 

= pk'<zk> 

- p.' 

It- 

,(Zk) 

= 0 

*k,3  (*kl- 

ak2‘  ak3-  ak4> 

= Pli'<zk) 

- pkl  l 

(Zk) 

= 0 

k=  1,  2 

’ » ■ • * 

. (n  - 

where  P£(x)  is  the  first  derivative  of  the  function  P^Cx) , with  respect 
to  the  coordinate  x,  and  P''(x)  is  the  second  derivative.  These  constraints 
express  the  conditions  required  for  the  ordinate,  slope,  and  second  deri- 
vatives of  the  fitting  polynomials  for  adjacent  intervals  to  be  continuous 
at  the  interval  boundary.  This  continuity  constrains  the  formulation  to 
provide  a more  nearly  correct  representation  of  the  real  physical  data. 


Lagrange's  method  of  undetermined  multipliers  (Refs.  12  and  13)  can 
be  used  to  minimize  Eq.  (6)  subject  to  the  constraints  expressed  by 
Eqs.  (7).  Let 
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3 

F,  (an.  ai2.  ai3,  ai4)  = SI+^  A], i^l,! 

j=l 


3 

Makl’  ak2’  ak3’  “k^  = Sk  (Ak-  l.j  &Y-  l.j  + Akj^kj^ 

j-1 


(8) 


|Hn  (an]’  an2>  an3’  an4^  = Sn  An- l.j  ^n-  1 ,j 


k = 2,  3,  . . . , (n  — 1 ) 


where  n is  the  number  of  intervals  and  the  A values  are  the  Lagrange 
multipliers. 


The  parameters,  A . , and  the  coefficients,  a 
mj 

necessary  minimization  are  determined  from  the  4n 


which  effect  the 
equations 


(«) 


k = 1,  2,  . . . , n 

i = l,  2,  3,  4 


(9) 


and  the  3(n-1)  equations 


W)*mJ  - 0 ; m=1'2>  ' • ■(n“1) 

m’J  j = 1,  2,  3 


(10) 


The  constant  multiplier,  1/2,  is  introduced  in  Eqs.  (9)  and  (10)  to  put 
the  equations  in  a form  which  is  convenient  for  subsequent  developments. 
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By  expanding  Eq.  (9)  for  k = 1,  one  obtains 


mi 

mJ 

mI 

an  ml 

* al 

i2  E*'  - 

i=  i 

a13 

Zi' 

i=  i 

a14  J]x<i  + 
i=  I 

'*  A1 1 

i=  i 

:n! 

mi 

ml 

2 

al>£ 

i=  1 

9 

“12 

i=  1 

■\  a 

13 

i=  1 

- auE'i 

1=  I 

/1 

2 

A 1 1 -*•  ’/.•  | A 

+ A13 

i=  1 


m] 


all 


1 


2Z'j  A 12 


(11) 


+ 6Z,  Aj3  = I; 

i=  1 


m, 


aJl 


i=  I 


2 


A 


II 


+ 3 Z 1 A]2  + 15Zj  A j 3 = 

i=  1 


where  the  x^’s  are  the  displacement  values,  is  the  second  endpoint, 
m1  is  the  number  of  points  in  the  first  interval,  and  Ii  = I(x.). 

The  constant  multiplier,  1/2,  in  Eq.  (9)  has  the  effect  of  eliminating 
the  factor  of  2,  which  is  introduced  in  Eq.  (9)  by  the  differentiation. 
Expanding  Eq.  (9)  for  k = 2,  3,  ....  n yields  more  equations  similar 
in  form  to  Eq.  (11). 


17 


AEDC-TR-76-163 


Expanding  Eq.  (10)  for  m = '1  yields 


'2all  + ~ al  2 


a13  + 


a14  — ^a-n  " a<>9  ~ 


*2 1 


22 


a2  3 ~ 


J24 


= 0 


3 5 

— 2Z  -i  a no  — 3Z  i at,  A — 0 


1 “ 1 2 + + ° \ “14  " ^ 1 “22  * ^]a23  ~ 0£jla24 

2 4 2 4 

a 1 2 + ^ ^ 1 a]  3 "t"  1 3 Z j s j ^ “ a22  — 6 Z 1 ajs  — 1 S Z j & 2 ^ = 0 


(12) 


Equations  of  similar  form  are  obtained  for  m = 2,  3,...,  (n-1). 


Equations  (4)  and  (10),  with  Eqs.  (11)  and  (12)  describing  represen- 
tative algebraic  details,  are  the  least-squares  spline  fit  equations  and 
represent  a system  of  equations  which  can  be  written  in  the  matrix- 
vector  notation 

BA  = C (13) 


where  B is  a (7n  - 3)  by  (7n  - 3)  symmetric  matrix,  A is  a (7n  - 3)  by 

1 matrix  which  has  as  elements  the  4n  coefficients,  a,  . , and  the  3n  - 3 

ki 

multipliers,  X , and  C is  a (7n  - 3)  by  1 matrix,  where  n is  the  number 
of  intervals  into  which  the  data  points  have  been  divided. 


The  matrix  B can  be  partitioned  into 


B = 


(14) 


where  R is  a 4n  x 4n  block  diagonal  matrix  which  is  defined  as 


R = 


(15) 
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The  matrices  P are  4 by  4 symmetric  matrices  which  are  defined  (letting, 


N is  a 4n  by  (3n  - 3)  block  band  matrix  defined  as 


(17) 
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The  matrix  N consists  of  a block  diagonal  and  one  block  lower  co-diagonal. 
T 

N is  the  transpose  of  the  matrix  N.  The  matrices  are  4 by  3 matrices 
defined  as 


\r 

'2 


z r 


Qj  = 


z4 

J 


(18) 


2 7.  6Z7 


Z6 

] 


3z: 


15Z4 


Noting  the  right  side  of  Eqs.  (11)  and  (12),  one  can  write  the  matrix 
C as 

C = GY  (19) 


where  G is  a (7n  - 3)  by  p matrix  of  certain  powers  of  x,  the  p-indepen- 
dent  data  points,  with  the  elements  of  the  last  3n  - 3 rows  all  zero,  and 
y is  a p by  1 matrix  of  the  p-dependent  variable  data  points  I(x)  (that  is, 
the  radiance  values).  Note  that  p = m^  + m2  + ...  + mn> 


For  ease  of  representation  the  matrix  G may  be  partitioned  into  a 
block  diagonal  matrix 


G = 


(20) 
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where  each  S^.  is  a 4 by  m^  matrix: 


mj- I*  1 


i-  1 


+ 1 


Writing  Eq.  (13)  as 


J- 


. + 2 

J- 


n'  i - 1 * 2 


A = B"  1 C 


1 

2 

Xmj-  1 + mj 

4 

‘"'j-  I •""j 
fi 

‘"'j-  , Mn j 


(21) 


(22) 


and  substituting  Eq.  (19)  yields 


A = B"  1 C Y 


(23) 


where  B is  a (7n  - 3)  by  (7n  - 3)  matrix,  which  is  the  inverse  of  B. 
Equation  (23)  can  be  written  as 


A = W Y 


(24) 


where  A is  the  (7n  - 3)  by  1 matrix,  which  has  as  elements  the  coefficients, 
aki’  an<^  ^ is  P by  1 matrix  of  data  intensity 


.-1 


values,  and  W equals  B G,  a (7n  - 3)  by  p matrix. 


Performing  the  integration  of  Eq.  (3)  with  the  polynomials  of  the 
form  of  Eq.  (4)  yields 


(r)  = - 


. 2 2 Ji 

tx  — r ) 


[ 2 * 2,j-  1 + J-  »3,j-  1 (*2+  2,2)  + -jV  “*■)-> 


2ar 


i = J 


(25) 


+ ~T~ a3 i (x2  + 2r2)  + «-  a4.  (3x4  + 4x2r2 
3 15 


♦ 8,4)] 


z. 


z. 


i-  ] 
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where  k is  chosen  so  that  r is  always  less  than  Z^.  Expanding  Eq.  (25) 
for  several  values  of  r,  one  can  write  the  resulting  system  of  equations 
in  the  matrix-vector  form 

E = MA  (26) 

where  E is  an  m by  1 matrix  of  values  of  the  emission  coefficient  e(r) 
evaluated  at  m different  values  of  r,  A is  a (7n  -3)  by  1 matrix  defined 
by  Eq.  (13),  and  M is  an  m by  (7n  - 3)  matrix  with  the  elements  of  the 
last  3n  - 3 columns  all  zero. 


Substituting  into  Eq.  (26),  the  expression  for  the  matrix  A from 
Eq.  (24)  yields 

E = MWY  (27) 

By  performing  the  matrix  multiplication 


F = M ff 


(28) 


Eq.  (22)  can  bfe  written  as 

E = FY  (29) 

where  E is  the  m by  1 matrix  of  values  of  the  emission  coefficient  e(r) 
evaluated  at  m different  values  of  r,  Y is  a p by  1 matrix  of  the 
p-dependent  variable  data  points  I(x),  and  F is  an  m by  p matrix  providing 
the  transformation  from  the  Y space  to  the  E space. 


2.2  ERROR  ANALYSIS  TECHNIQUE 


To  complete  the  development  of  the  technique  to  include  the  error 
propagation  analysis,  it  is  necessary  to  utilize  some  definitions  and 
results  from  statistics  theory  (Refs.  14  and  15).  The  expected  value 
of  a random  variable  x,  denoted  by  6 (x) , is  obtained  by  finding  the 
average  value  of  the  function  over  all  possible  values  of  the  variable. 
The  expected  value  of  a matrix  or  vector  M,  denoted  by  £(M) , is  the 
matrix  of  the  expected  values  of  the  elements  of  M.  The  moments  of  a 
distribution  are  the  expected  values  of  the  powers  of  the  variable  which 
have  the  given  distribution. 
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Let  yE  denote  the  first  moment  or  mean  of  the  emission  coefficient 
function  e(r)  and  let  denote  the  mean  of  the  radiance  function  I(x). 
Then 

Ht  = S(F.)  = fi(KY)  = V&{\)  = l’>,  (30) 

The  covariance  matrix  of  E,  denoted  by  [E]  , is  the  m by  m 

symmetric  matrix  defined  as 

IE]C,  = & [(F.  - /if)T]  (31) 

where  m is  the  number  of  different  values  of  r at  which  the  emission 
coefficient  e(r)  is  evaluated.  The  following  steps  yield  an  expression 
which  can  be  used  to  evaluate  the  covariance  matrix  of  E: 

fKlcv  = 6 [<E-|xf)(E 

= S [(FY  - F/i,)(  FY  - F/t,)1] 

= & [F(  Y -/£j)(F(Y  -^,))T]  (32) 

= 6 [f(Y  -Mi)(Y  -M[)'r  ft] 

= I*'  [(  Y - fij)  ( Y - f4, )TJ  ) F 1 

[EJCV  = F[Y]„FT 

The  term  [Y]  in  Eq.  (32)  denotes  the  covariance  matrix  of  Y,  which  is 
an  m by  m symmetric  matrix. 

The  diagonal  elements  of  the  matrix  fY J are  identified  with  the 
squares  of  the  standard  deviations  of  the  radiance  measurements.  The 
nondiagonal  elements  are  the  covariances  between  the  various  elements  of 
the  matrix  Y.  In  the  present  problem,  since  the  individual  observations 
are  independent,  then  the  elements  of  Y (that  is,  the  radiance  data) 
are  assumed  to  be  uncorrelated.  Therefore,  all  of  the  nondiagonal 
elements  of  the  matrix  [Y] ^ are  zero,  and  thus  [Y]  is  simply  a diagonal 
matrix  containing  the  variances  of  the  radiance  data. 
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The  matrix  F provides  the  transformation  from  the  Y space  to  the 
E space.  The  form  of  the  matrix  is  determined  by  the  raw  data  and  the 
choice  of  numerical  technique.  The  [e1cv  is  the  covariance  matrix  of  E, 
and  the  diagonal  elements  are  the  variances  of  the  emission  coefficient 
values  which  have  been  calculated  by  the  least-squares  spline  fit  approxi- 
mation to  the  data.  The  off-diagonal  elements  of  the  matrix  [E]  give  the 
covariance  between  the  various  emission  coefficient  values.  The  standard 
deviations  of  the  emission  coefficient  values  are  the  positive  square 
roots  of  the  diagonal  elements  of  [e]^.  Thus,  determination  of  the 
matrix  F,  Eq.  (28),  and  its  use  in  Eqs.  (29)  and  (32)  provides  the  formal 
solution  to  the  inversion  problem  and  the  propagation  of  the  random 
errors  associated  with  the  measurements. 

3.0  RESULTS  AND  DISCUSSION 


3.1  INTRODUCTION 

In  the  determination  of  the  solution  to  a specific  inversion  pro- 
blem, described  formally  by  Eqs.  (27)  and  (30),  there  are  generally  four 
parameters  which  may  be  varied  by  the  user.  These  four  parameters  are 
(1)  the  total  number  of  data  p'oints,  (2)  the  distribution  of  the  data 
points  chosen  for  analysis,  (3)  the  number  of  intervals  into  which  the 
data  are  divided,  aijd  (4)  the  number  of  points  in  each  interval.  Variation 
of  each  of  these  parameters  affects  the  elements  of  the  matrix  F,  Eq.  (28), 
and  the  subsequent  inversion  and  random  error  propagation  results.  Included 
in  this  section  are  results  of  the  Abel  inversion  of  several  sets  of 
analytic  data  chosen  to  illustrate  the  effect  of  variations  of  the  above 
four  parameters  on* the  results  of  the  inversion.  Also  included,  as  an 
illustration  of  the  application  of  the  method  to  typical  experimental 
data,  are  the  results  of  the  inversion  of  data  taken  in  a recent  research 
experiment. 
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3.2  ANALYTIC  CONTINUOUS  TEST  DATA 

Thirty-one  data  points  are  used  to  illustrate  the  technique.  These 
data  points  were  generated  from  the  function 

I(x)  = e"x  ; 0<x<3.0  (33) 

which  may  be  inverted  by  Eq.  (2)  to  yield  the  analytic  function 

f { r)  = JL  e“r  erf  [(R2-  r2)S]  (34) 

yjn 

The  number  of  data  points  and  the  displacement  distribution  were  fixed, 
leaving  only  the  number  of  intervals  and  the  number  of  points  per  inter- 
val to  be  varied.  The  inversion  was  applied  to  several  cases  with 
different  values  of  the  free  parameters.  A standard  deviation  of  10 
percent  of  each  radiance  value  was  assumed  in  each  case.  The  displace- 
ment radiance  and  standard  deviations  used  are  listed  in  Table  1. 

The  initial  data  configuration  examined  was  formed  by  dividing  the 
data  into  four  intervals,  with  seven  points  in  the  first  interval  and 
eight  points  in  the  remaining  three  intervals.  The  results  of  inverting 
the  data  with  this  configuration  of  the  test  data  are  presented  in 
Table  2 and  Figs.  2 and  3.  Figure  2 displays  the  input  radiance  data 
as  points  and  the  results  of  the  curve  fit  as  a continuous  line.  Error 
intervals  equal  in  magnitude  to  the  radiance  standard  deviations  are 
shown  for  each  data  point.  Figure  3 displays  the  profile  of  the  calculated 
emission  coefficients  along  with  the  calculated  error  interval  for  each 
emission  coefficient  value. 

Another  configuration  of  the  same  set  of  test  data  was  inverted 
after  distributing  the  31  data  points  into  six  intervals  with  five  points 
in  the  first  five  intervals  and  six  points  in  the  sixth  interval.  The 
results  of  this  case  are  presented  in  Table  3 and  Figs.  4 and  5.  Figure 
4 shows  the  radiance  and  Fig.  5 the  emission  coefficient  profile  for 
this  configuration  of  the  data.  , 
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A comparison  of  Tables  1 and  2 reveals  that  the  percentage  errors 
between  the  curve  fit  data  and  the  input  data  are  generally  smaller  in 
magnitude  for  the  six-interval  configuration  than  for  the  four- interval 
configuration.  Of  the  31  percentage  errors,  24  are  smaller  for  the 
six- interval  configuration.  This  suggests  that  the  curve  fit  of  the 
radiance  data  was  better  for  the  case  of  six  intervals.  Furthermore,  of 
the  31  standard  deviation  intervals  displayed  in  each  of  Figs.  3 and  5, 

22  are  smaller  for  the  six- interval  case.  This  fact  is  illustrated  more 
clearly  by  comparing  the  calculated  standard  deviations  of  the  emission 
coefficient  values  for  the  two  cases  which  are  listed  in  the  last  columns 
of  Tables  2 and  3. 

Table  1.  Checkout  Data  with  e*2  Values 
Used  as  Input  Data 


Displacement 

Radiance 

(Data) 

Standard  Deviation 

0.0 

1.0000 

E-00 

1.0000 

E-01 

0.  1 

D.  0005 

E-01 

9.  BOOS 

E-02 

0.2 

B.  607  D 

E-01 

0.  6079 

E-02 

0.3 

0.  1393 

E-01 

9. 1393 

E-02 

0.4 

8.5214 

E-01 

8.  5214 

E-02 

0.5 

7.7880 

E-01 

7. 7BB0 

E-02 

0.6 

G.97GB 

E-01 

G. 976B 

E-02 

0.7 

6.  1263 

E-01 

C. 1263 

E-02 

0.8 

5. 2726 

E-0- 

5.2720 

E-02 

o.  o 

4.4863 

E-01 

4.4860 

E-02 

1.0 

3.6788 

E-01 

3.6788 

E-02 

1. 1 

2.  9820 

E-01 

2.9820 

E-02 

1.2 

2.3693 

E-01 

2.3693 

E-02 

1.3 

1.8452 

E-01 

1.8452 

E-02 

1.4 

1.4086 

E-01 

1.4086 

E-02 

1.5 

1.0540 

E-01 

1.  0540 

E-02 

1.6 

7. 7305 

E-02 

7.7305 

E-03 

1.7 

5. 5576 

E-02 

5.5576 

E-03 

l.B 

3.9164 

E-02 

3.9164 

E-03 

1.9 

2. 7052 

E-02 

2.  7052 

E-03 

2.0 

1.8316 

E-02 

1.8316 

E-03 

2.  1 

1.2155 

E-02 

1.2155 

E*-03 

2.2 

7.0C71 

E-0; 

7.  BC71 

E-04 

2.3 

5.0418 

E-05 

5.0418 

E-04 

2.4 

3. 1511 

E-0  3 

3.1511 

E-04 

2.5 

1.0305 

E-03 

1.9305 

E-04 

2.6 

1.  1592 

E-03 

1.  1592 

E-04 

2.7 

6.8233 

E-04 

6.  8233 

E-05 

2.8 

3.93G7 

E -04 

3.0367 

E-05 

2.9 

2.2263 

E-04 

2.2263 

E-05 

3.0 

1.2341 

E-04 

1,2341 

E-05 
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Table  2. 


Inversion  Results  Using  Four  Intervals  of 
e*2  Values  as  Input  Data 


Displacement 

Radiance 

(Calculated) 

Percent  Error 
between  Radiance 
Data  and 

Calculated  Radiance 

Emission  Coefficient 

Standard  Deviation 
(Emission  Coefficient) 

0.0 

1.000687 

E-00 

6.873  = 65 

E-02 

5. 675450 

E-01 

1.  277898 

E-01 

1.000000 

E-01 

9. 904464 

E-01 

4.003959 

E-02 

5.613535 

E-01 

1. 183743 

E-01 

2.  000000 

E-01 

9.  604699 

E-01 

-3.331861 

E-02 

5. 432904 

E-01 

9.255000 

E-02 

3.000000 

E-01 

9. 128938 

E-01 

-1. 133770 

E-Ol 

5.  148417 

E-01 

5. 781569 

E-02 

4.000000 

E-01 

8.509365 

E-01 

-1.413484 

E-01 

4. 783169 

E-01 

3.000306 

E-02 

5. 000000 

E-01 

7.  783793 

E-01 

-5.401784 

E-02 

4.365562 

E-01 

3.  155638 

E-02 

G.  000000 

E-01 

6.  968567 

E-01 

1.686561 

E-01 

3.923847 

E-01 

3. 612286 

E-02 

7. 000000 

E-01 

6. 149943 

E-Ol 

3.859220 

E-01 

3.470961 

E-01 

2.810619 

E-02 

8. 000000 

E-01 

5.293764 

E-01 

3.956902 

E-01 

3.001493 

E-01 

1.929250 

E-02 

9. 000000 

E-01 

4.455566 

E-01 

-6. 784139 

E-01 

2.  531834 

E-01 

1,210634 

E-02 

1.  000000 

E-01 

3.668708 

E-00 

-2.743241 

E-01 

2.081161 

E-01 

8.488109 

E-03 

1.  100000 

E-00 

2. 960949 

E-00 

-7. 059348 

E-01 

1.  668016 

E-01 

9.052147 

E-03 

1.200000 

E-00 

2. 351273 

E-00 

-7. 608546 

E-01 

1.308842 

E-01 

1.032227 

E-03 

1.  300000 

E-00 

1.846380 

E-00 

6.  396129 

E-02 

1.015773 

E-01 

9.  970046 

E-03 

1.400000 

E-00 

1.436843 

E-01 

2.005025 

E-00 

7. 925596 

E-02 

7. 552663 

E-03 

1. 500000 

E-00 

1.094268 

E-01 

3.820451 

E-00 

6.210446 

E-02 

4. 133476 

E-03 

1.  600000 

E-OO 

8.000591 

E-02 

3.493838 

E-00 

4.684850 

E-02 

2.054649 

E-03 

1.  700000 

E-00 

5.579709 

E-02 

3. 978129 

E-01 

3.352474 

E-02 

2.698357 

D-03 

1.800000 

E-00 

3.  703233 

E-02 

-5.442925 

E-00 

2. 244468 

E-02 

3.813321 

E-03 

1.900000 

E-00 

2.362409 

E-02 

-1.267157 

E-01 

1.  384293 

E-02 

4.234182 

E-03 

2.000000 

E-00 

1.  510630 

E-02 

-1.752401 

E-01 

7.  840132 

E-03 

3,858621 

E-03 

2.  100000 

E-00 

1.057220 

E-02 

-1.057220 

E-01 

4.382294 

E-03 

2. 596808 

E-03 

2. 200000 

E-00 

8.606659 

E-03 

8.872517 

E-00 

3. 113924 

E-03 

8.552391 

E-04 

2.  300000 

E-00 

7.346027 

E-03 

4.570248 

E-01 

2. 879137 

E-03 

9. 776775 

E-04 

2.400000 

E-00 

5. 995310 

E-03 

9.026087 

E-01 

2. 571927 

E-03 

1.831465 

E-03 

2. 500000 

E-00 

4.613937 

E-03 

1.390022 

E-02 

2.  192267 

E-03 

1.821228 

E-03 

2.  600000 

E-00 

3.264921 

E-03 

1.816530 

E-02 

1.  750231 

E-03 

1.024300 

E-03 

2. 700000 

E-00 

2.026027 

E-03 

1.969277 

E-02 

1.  263463 

E-03 

3.  672080 

E-04 

2.  800000 

E-00 

9.912250 

E-04 

1.517908 

E-02 

7.616202 

E-04 

2.026932 

E-03 

2.  900000 

E-00 

2.722152 

E-04 

2.227249 

E-01 

2. 972067 

E-04 

3.  192960 

E-03 

3.000000 

E-00 

0.00 

-1.000000 

E-02 

0.0 

0.0 

Number  of  points:  31 

Number  of  Intervals:  4 

Number  of  points  per  interval:  7 8 8 8 
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Number  of  Points:  31 

Number  of  Intervals:  4 

Number  of  Points  per  I ntervai:  7 8 8 8 

$ Data  and  Uncertainty 


0 0.60  1.20  1.80  2.40  3.00 


Displacement,  arbitrary  units 


Figure  2.  Intensity  for  four-interval 


Radius,  arbitrary  units 


Figure  3.  Emission  coefficient  profile  for  four-interval 
e-*2  test  data. 
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Table  3.  Inversion  Result  with  Six-Interval 
e-*2  Values  Used  as  Input  Data 


Displacement 

Radiance 

(Calculated) 

Percent  Error 
tetween  Radiance 
Data  and 

Calculated  Radiance 

Emission  Coeffic.ent 

Standard  Deviation 
(Emission  Coefficient) 

0.0 

1.000057 

E-00 

5.  727384 

E-03 

5. 644821 

E-01 

1.887621 

E-01 

1.  000000 

E-01 

9.90060B 

E-01 

1.095800 

E-03 

5.  587444 

E-01 

1.589008 

E-01 

2.000000 

E-01 

9.607041 

E-01 

-8.935534 

E-03 

5.419311 

E-01 

8. 748871 

E-02 

3. 000000 

E-01 

9. 138023 

E-01 

-1.3977G1 

E-02 

5. 151949 

E-01 

4. 705685 

E-02 

4.000000 

E-01 

8.521048 

E-0l' 

-4. 125894 

E-03 

4.802804 

E-01 

6.393484 

E-02 

5. 000000 

E-01 

7.788967 

E-01 

1.241464 

E-02 

4.390734 

E-01 

4.259269 

E-02 

6,000000 

E-01 

6. 978503 

E-01 

2.441074 

E-02 

3. 932836 

E-01 

2.312669 

E-02 

7.  000000 

E-01 

6. 129601 

E-01 

5. 388051 

E-02 

3. 452108 

E-01 

2.  111331 

E-02 

B.  000000 

E-01 

5. 260066 

E-01 

1. 359116 

E-01 

2.  972415 

e-o: 

2.471700 

E-02 

9. 000000 

E-01 

4.460712 

E-01 

-5. 636391 

E-01 

2. 514763 

E-01 

2.04G078 

E-02 

1.000000 

E-00 

3.690464 

E-01 

3.  1705 82 

E-01 

2.088244 

E-01 

1.096933 

E-02 

1. 100000 

E-00 

2.  935861 

E-01 

1.294751 

E-01 

1.  690072 

E-01 

7. 882454 

E-03 

1„  200000 

E-00 

2.366441 

E-01 

-1.20iM05 

E-01 

1.335991 

e-o: 

9. 789215 

E-03 

1.300000 

E-00 

1.842825 

E-01 

-1.287393 

E-01 

1.035031 

E-01 

0.862244 

E-03 

1.400000 

E-00 

1.413117 

E-01 

3.206904 

E-01 

7.938548 

E-02 

6. 747721 

E-03 

1. 500000 

E-00 

1.060471 

E-01 

6. 139940 

E-01 

6.027553 

E-02 

2.732949 

E-03 

1. 600000 

E-00 

7.  724082 

E-02 

-8.302233 

E-02 

4.413603 

E-02 

2.944425 

E-03 

1.  700000 

E-00 

5.489177 

E-02 

-1. 231169 

E-00 

3. 111812 

E-02 

2.669870 

E-03 

1,800000 

E-00 

3. 854638 

E-02 

-1.577005 

E-00 

2. 141045 

E-02 

3.67C396 

E-03 

1. 500000 

E-00 

2.711363 

E-02 

2.278071 

E-01 

1.483089 

E-02 

1.  931244 

E-03 

2.000000 

E-00 

1.881522 

E-02 

2. 725-04 

E-00 

1. 060456 

E-02 

6.  662148 

E-04 

2.  100000 

E-00 

1.246579 

E-02 

2.  556850 

E-00 

7.  207111 

E-03 

1. 797132 

E-03 

2.200000 

E-00 

7.899407 

E-03 

-9. 729131 

E-02 

4. 604004 

E-03 

8.499760 

E-04 

2.300000 

E-00 

4.886028 

E-03 

-3.089606 

E-00 

2.  794489 

E-03 

9.  729830 

E-04 

2.400000 

E-00 

3.061859 

E-03 

-2.832055 

E-00 

1. 712425 

E-03 

4.271174 

E-04 

2.500000 

E-00 

1.917642 

E-03 

-6.660687 

E-01 

1.  137634 

E-03 

7.384108 

E-04 

2. 600000 

E-00 

1.084741 

E-03 

-6.423344 

E-Gu 

7.08G60B 

E-04 

1.008421 

E-03 

2.  700000 

E-00 

5.  185475 

E-04 

-2.400342 

E-Oi 

3.841813 

E-04 

3.792780 

E-04 

2.  600000 

E-00 

1.  843718 

E-04 

-5.  316590 

e-o; 

1.  633415 

E-04 

3.629174 

E-04 

2.  900000 

E-00 

3.313120 

E-05 

-8.510728 

E-01 

3. 9932 10 

E-O  5 

1.  336839 

E-03 

3.000000 

E-00 

0.0 

-l.dooooo 

E-02 

0.0 

0.0 

Number  of  points.  31 

Number  of  intervals:  6 

Number  of  points  per  interval:  5 5 5 5 5 6 
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Figure  4.  Intensity  for  six-interval  e*2  test  data. 


Figure  5.  Emission  coefficient  profile  for  six-interval 
e*2  test  data. 
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The  results  of  this  comparison  point  to  the  conclusion  that  the 
emission  coefficient  values  obtained  using  the  six-interval  configuration 
are  more  accurate  for  the  particular  set  of  data  tested  than  are  the 
values  obtained  using  the  four-interval  configuration.  A comparison  of 
the  emission  coefficient  values  obtained  for  the  two  cases  with  the 
analytically  determined  values  obtained  by -evaluating  Eq.  (1)  at  the 
proper  displacement  values  shows  that  this  conclusion  is  indeed  correct. 
Table  4 lists  the  analytical  values  of  the  emission  coefficients  for  this 
particular  problem  and  the  percent  error  between  the  analytic  emission 
coefficient  values  and  the  calculated  values  for  the  two  cases.  Of 
the  31  percentage  errors  listed  for  each  configuration,  26  are  smaller 
in  magnitude  for  the  six- interval  case. 

Table  5 and  Figures  6 and  7 show  the  results  of  inverting  the  data 
with  a configuration  of  six  intervals  with  an  uneven  dispersion  of  the 
data  points  among  the  intervals — two  points  in  the  first,  third,  and 
fourth  intervals,  eight  points  in  the  second  and  fifth  intervals,  and 
nine  points  in  the  sixth  interval.  A comparison  of  these  results  with  the 
results  of  the  six-interval  configuration  with  an  even  dispersion  of  the 
data  (Table  3,  Figs.  4 and  5)  suggests  that  the  configuration  of  data 
with  the  data  points  dispersed  uniformly  among  the  intervals  yields  more 
accurate  emission  coefficient  values  for  these  analytical  data.  This 
conclusion  is  supported  by  the  data  listed  in  Table  6,  which  shows  the 
percentage  error  between  the  analytic  emission  coefficient  values  and  the 
calculated  values  for  the  two  cases.  Examination  of  the  standard  deviations 
of  the  emission  coefficients  for  the  six-interval  configuration  with 
an  even  and  uneven  distribution  of  the  data  points  among  the  intervals 
reveals  that  the  even  distribution  provides  noticeably  smaller  standard 
deviations  near  the  centerline.  Standard  deviations  are  also  generally 
smaller  for  larger  displacement  values  for  the  even  distribution,  although 
the  differences  are  not  as  great. 
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Table  4.  Percentage  Error  Between  Analytic  Emission 
Coefficient  Values  and  Calculated  Values 


Displacement 

Analytic  Emission 
Coefficient 

Percent  Error 
(Four  Intervals) 

Percent  Error 
(Six  Intervals) 

0.0 

5.  634561 

E-01 

7.256821 

E-01 

1.820905 

E-01 

0.1 

5.  578497 

E-01 

6.280903 

E-01 

1.603831 

E-01 

0.2 

5.413085 

E-01 

3.  661313 

E-01 

1.  150176 

E-01 

0.3 

5.  149086 

E-01 

-1.299259 

E-02 

5.560210 

E-02 

0.4 

4,800495 

E-01 

-3.609211 

E-01 

4.809920 

E-02 

0.5 

4.387322 

E-01 

-4. 959745 

E-01 

7.776954 

E-02 

0.  6 

3.  929919 

E-01 

-1.545070 

E-01 

7.425089 

E-02 

0.  7 

3.  450153 

E-01 

-6.031037 

E-01 

5.666415 

E-02 

0.8 

2.  969277 

E-01 

1.084978 

E-00 

1.056823 

E-01 

0.9 

2.  504572 

E-01 

1.088489 

E-00 

4.068959 

E-01 

1.0 

2.070764 

E-01 

5.020852 

E-01 

8.441329 

E-01 

1.1 

1.678024 

E-01 

-5. 964158 

E-01 

7.  716219 

E-01 

1.2 

1.  332712 

E-01 

-1,  791085 

E-00 

2.460397 

E-01 

1.3 

1.037396 

E-01 

-2.084354 

E-00 

-2.279747 

E-01 

1.4 

7.  915297 

E-02 

1.301404 

E-01 

2.937477 

E-01 

1.  5 

5.918566 

E-02 

4.  931634 

E-00 

1.841443 

E-00 

1.6 

4.  337465 

E-02 

8.008941 

E-00 

1.  755357 

E-00 

1.7 

3. 114230 

E-02 

7.  650174 

E-00 

-7.  764359 

E-02 

1.8 

2.  191468 

E-02 

2.418470 

E-00 

-2.300878 

E-00 

1.9 

1.  510670 

E-02 

-8.  365565 

E-00 

-1.494702 

E-00 

2.0 

1.020432 

E-02 

-2.316850 

E-01 

3. 922260 

E-00 

2.  1 

6.  746728 

E-03 

-3.504564 

E-01 

6.823827 

E-00 

2.2 

4.  368733 

E-03 

-2.872249 

E-01 

5.385360 

E-00 

2.3 

2.768276 

E-03 

4.004695 

E-00 

9.468709 

E-01 

2.4 

1.  714000 

E-03 

5.005408 

E-01 

-9.200700 

E-02 

2.5 

1.  103632 

E-03 

9.864114 

E-01 

3.080837 

E-00 

2,6 

6.  103361 

E-04 

1.  867651 

E-02 

1.610981 

E-01 

2.7 

3.  483524 

E-04 

2. 626968 

E-02 

1.028519 

E-01 

2.8 

1.  909872 

E-04 

2.  987808 

E-02 

-1.  447511 

E-01 

2.9 

9.  789680 

E-05 

2.035918 

E-02 

-5.917934 

E-01 

3.0 

3.481320 

E-05 

-1.000000 

E-02 

-1.000000 

E-02 
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Displacement 


Table  5.  Inversion  Results  with  Unevenly  Dispersed 
e*2  Values  Used  as  Input  Data 


Radiance 

(Calculated) 


Percent  Error 
between  Radiance 
Data  and 

Calculated  Radiance 


Emission  Coefficient 


Standard  Deviation 
(Emission  Coefficient) 


0.  00 

9.  999292 

E-01 

-7. 174047 

E-03 

5.  609455 

E-01 

8. 117026 

E-01 

1. 000000 

E-01 

9.  904682 

E-01 

4.224036 

E-02 

5.  588918 

E-0! 

1.611006 

E-01 

2. 000000 

E-01 

9.  610668 

E-01 

2.881206 

E-02 

5.431304 

E-01 

9.991750 

E-02 

3.000000 

E-01 

9.  138335 

E-01 

-1.055980 

E-02 

5.  161061 

E-01 

7.471066 

E-02 

4.000000 

E-01 

8. 516753 

E-01 

-5.453786 

E-02 

4.804729 

E-01 

4. 784840 

E-02 

3.000000 

E-01 

7. 781568 

e-o: 

-8.258385 

E-02 

4.382930 

E-01 

3. 147328 

E-02 

6. 000000 

E-01 

G.  972846 

E-01 

-5. 656698 

E-02 

3.920035 

E-01 

3.479315 

£-02 

7. 000000 

E-01 

6.  130743 

E-01 

7. 251822 

E-02 

3.442192 

E-01 

3.938439 

E-02 

8. 000000 

E-01 

5.290217 

E-01 

3.284228 

E-01 

2.974210 

E-01 

3.213198 

E-02 

9. 000000 

E-01 

4.474787 

E-01 

-2.499491 

E-01 

2.  533768 

E-01 

2.010979 

E-02 

1.000000 

E-00 

3. G90S51 

E-01 

3.275864 

E-01 

2. 112019 

E-01 

2.918001 

E-02 

1.  100000 

E-00 

2. 960012 

E-01 

-7.373671 

E-01 

1.683103 

E-01 

2. 209935 

E-02 

1.200000 

E-00 

2.339679 

E-01 

-1.250211 

E-00 

1.294129 

E-01 

8.216915 

E-03 

1.  300000 

E-00 

1.850101 

E-01 

2. 656128 

E-01 

1.011093 

E-01 

1.  103693 

E-02 

1.400000 

E-00 

1.443999 

E-01 

2. 513067 

E-00 

8,076948 

E-02 

7.790830 

E-03 

1.500000 

E-00 

1.  088006 

E-01 

3.  22 63 34 

E-00 

6.254847 

E-02 

4.444626 

£-03 

1.600000 

E-00 

7.871845 

E-02 

1. 628408 

E-00 

4. 633172 

E-02 

2.  182256 

E-03 

1.700000 

E-00 

5.457954 

E-02 

-1. 792976 

E-00 

3.252025 

E-02 

1.997622 

E-03 

1.  800000 

E-00 

3. 646029 

E-02 

-6. 903555 

E-00 

2. 143476 

E-02 

2.740380 

E-03 

1. 5)00000 

E-00 

2.40170G 

E-02 

-1.  121892 

E-01 

1.327453 

E-02 

2. 979927 

E-03 

2.  000000 

E-00 

1,641107 

E-02 

-1.040037 

E-01 

8,051507 

E-03 

2,499504 

E-03 

2.  100000 

E-00 

1.222921 

E-02 

6. 105477 

E-01 

5.459507 

E-03 

1,413835 

E-03 

2.  200000 

E-00 

9. 487126 

E-03 

1.  998237 

E-01 

4.376879 

E-03 

3.050172 

E-04 

2.  300000 

E-00 

7. 102601 

E-03 

4.087430 

E-01 

3.451966 

E-03 

7.71G312 

E-04 

2.  400000 

E-00 

5.054098 

E-03 

6.070390 

E-01 

2.  C 17529 

E-03 

1.  100082 

E-03 

2.  500000 

E-00 

3.383133 

E-03 

7.  524648 

E-01 

1.883582 

E-03 

9.  757774 

E-04 

2.  600000 

E-00 

2.061617 

E-03 

7. 784823 

E-01 

1.260428 

E-03 

4. 451294 

E-04 

2.  700000 

E-00 

1.090593 

E-03 

5. 983366 

E-01 

7.543510 

E-04 

4.054118 

E-04 

2.  800000 

E-00 

4.489278 

E-04 

1.403G59 

E-01 

3. 714054 

E-04 

1.  337966 

E-03 

2.  900000 

E-00 

1.016430 

E-04 

-5.420660 

E-01 

1. 149017 

:-:-04 

1.95  5461 

E-03 

3.  000000 

E-00 

0.  0 

-1.000000 

E-02 

0.0 

o.’o 

Number  of  points:  31 
Xumoer  of  intervals:  6 

Number  of  points  per  interval:  2 8 2 2 8 9 
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Figure  6.  Intensity  for  unevenly  dispersed  e*2 
test  data. 


Radius,  arbitrary  units 

Figure  7.  Emission  coefficient  profile  for  unevenly 
dispersed  e'*2  test  data. 
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Table  6.  Percentage  Error  Between  Analytic  Emission 
Coefficient  Values  and  Calculated  Values, 
Even  and  Uneven  Data 


Displacement 

Percent  Error 
(Even  Dispersion) 

Percent  Error 
(Uneven  Dispersion) 

0.0 

1.  820905 

E-01 

-4.455715 

E-01 

0. 1 

1.603837 

E-01 

1.858055 

E-01 

0.2 

1. 150176 

E-01 

3.365733 

E-01 

0.  3 

5.560210 

E-02 

2.325655 

E-01 

0.4 

4.809920 

E-02 

8.819924 

E-02 

0.5 

7. 776954 

E-02 

-1.001066 

E-01 

0.6 

7.425089 

E-02 

-2.515065 

E-01 

0.7 

5.  666415 

E-02 

-2.307434 

E-01 

0.8 

1.056823 

E-01 

-1.661347 

E-01 

0.9 

4.068959 

E-01 

1. 165708 

E-00 

1.0 

8.441329 

E-01 

-1.992260 

E-00 

1.1 

7.716219 

E-01 

3.026774 

E-01 

1.2 

2.460397 

E-01 

-2.895074 

E-00 

1.3 

-2.279747 

E-01 

-2.535483 

E-00 

1.4 

2. 937477 

E-01 

2.042261 

E-00 

1.5 

1.841443 

E-00 

5.681799 

E-00 

1.6 

1. 755357 

E-00 

6.817507 

E-00 

1.7 

-7. 764359 

E-02 

4.424689 

E-00 

1.8 

-2.300878 

E-00 

-2. 189948 

E-00 

1.9 

-1.494702 

E-00 

-1.212814 

E-01 

2.0 

3. 922260 

E-00 

-2.109521 

E-01 

2.1 

6.823827 

E-00 

-1.907919 

E-01 

2.2 

5.385360 

E-00 

1.864614 

E-01 

2.3 

9.468709 

E-01 

2.469732 

E-01 

2.4 

-9. 200700 

E-02 

5.271464 

E-01 

2.5 

3.080837 

E-00 

7.070746 

E-01 

2.6 

1.610981 

E-01 

1.065138 

E-02 

2.7 

1.028519 

E-01 

1.165482 

E-02 

2.8 

-1.447511 

E-01 

9. 446612 

E-01 

2.9 

-5.917934 

E-01 

1.839171 

E-01 

3.0 

-1.000000 

E-02 

-1.000000 

E-02 
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As  a further  examination  of  the  technique,  the  31  ordinates  generated 
from  the  function  of  Eq.  (33)  were  randomly  scattered  by  a fractional 
amount  bounded  by  ±10  percent  of  the  radiance  values.  The  scattered  data, 
listed  in  Table  7,  were  tested  for  all  the  same  parameter  configurations 
used  for  the  unscattered  data.  Table  8 lists  the  results  of  inverting  the 
scattered  data  in  the  six-interval  configuration,  and  Figs.  8 and  9 display 
the  results  graphically.  A comparison  of  this  data  with  that  shown  in 

Table  7.  Checkout  Data  with  Scattered  e*2 
Values  Used  as  Input  Data 


Dia  placement 

Intensity  (Data) 

Standard  Deviation 

0.0 

9.  90000 

E-01 

9.  90000 

E-02 

0.  1 

9.  30650 

E-01 

9.  30650 

E-02 

0.2 

1.01844 

E-00 

1.01844 

E-01 

0.3 

9.  23070 

E-01 

9.  23070 

E-02 

0.4 

B. 60660 

E-01 

8.  60660 

E-02 

0.5 

7.24284 

E-01 

7.24284 

E-02 

0.6 

6.  97680 

E-01 

6.97680 

E-02 

0.7 

5.  88120 

E-01 

5.88120 

E-02 

0.8 

5.43110 

E-01 

5.43110 

E-02 

0.9 

4. 08230 

E-01 

4.08230 

E-02 

1.0 

3.45810 

E-01 

3.45810 

E-02 

1. 1 

2.  98200 

E-01 

2. 98200 

E-02 

1.2 

2.  55880 

E-01 

2.55880 

E-02 

1.3 

1.  95590 

E-01 

1.95590 

E-02 

1.4 

1.31000 

E-01 

1.31000 

E-02 

1.5 

1.  10670 

E-01 

1.  10670 

E-02 

1.6 

7.49860 

E-02 

7.49860 

E-03 

1.7 

5.89110 

E-02 

5.89110 

E-03 

1.8 

3.56390 

E-02 

3.  56390 

E-03 

1.9 

2.  62400 

E-02 

2.  62400 

E-03 

2.0 

1.  74000 

E-02 

1.  74000 

E-03 

2.  1 

1.  13040 

E-02 

1.  13040 

E-03 

2.2 

8.06520 

E-03 

8. 06520 

E-04 

2.3 

4.89050 

E-03 

4. 89050 

E-04 

2.4 

2.86750 

E-03 

2.  86750 

E-04 

2.5 

1.81470 

E-03 

1.81470 

E-04 

2.6 

1. 10120 

E-03 

1. 10120 

E-04 

2.7 

7.30090 

E-04 

7.  30090 

E-05 

2.8 

3.93670 

E-04 

3.93670 

E-05 

2.9 

2.  20400 

E-04 

2.20400 

E-OS 

3.0 

1.25880 

E-04 

1.25880 

E-05 
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Table  8.  Inversion  Results  with  Randomly  Scattered 
e-x^  Values  Used  as  Input 


Percent  brror 


Disidacument 

Radiance 

(Calculated) 

between  Radiance 
Data  and 

Calculated  Radiance 

Emission  Coefficient 

Standard  Deviation 
(Emission  Coefficient) 

C.O 

9.  739G92 

E-01 

-1.  114224 

E-00 

4. 84 1228 

K-01 

2.205469 

E-Cl 

1.000000 

Y -01 

9. 761542 

E-01 

4.  B8D508 

L-00 

4. 952502 

E-01 

1.  904222 

E-01 

2.000000 

K-01 

9.  627751 

E-01 

-5,465700 

E-00 

5.203799 

i-r-oi 

1.042947 

E-01 

3.000000 

B-0I 

9.  271917 

E-01 

4.465221 

E-01 

5.371826 

E-01 

5. 780840 

E-02 

4.000000 

E-01 

8. 605448 

E-01 

-1.338141 

E-02 

5.  181171 

L-Cl 

7. 889482 

K-02 

5. 000000 

E-01 

7,715935 

E-01 

5.  531891 

E-30 

4.59.635 

E-Cl 

5,243160 

E-02 

G.OOOOOO 

E-01 

G. 781933 

E-0*1 

-2.73307G 

L-00 

3.938677 

E-0 1 

2.  803216 

E-02 

7.  000000 

E-01 

5.  875534 

E-01 

-9.G33717 

E-02 

J. .114019 

E-01 

2.527556 

E-02 

R. 000000 

E-01 

5.044674 

E-01 

-7.  1 15067 

E-00 

2. 764865 

E-01 

2.979423 

K-02 

9. 000000 

E-01 

4.  JO '>600 

E-0) 

5.4G995B 

E-00 

2.  327.(88 

E-01 

2.474729 

E-02 

1,000000 

E-30 

3.  628808 

e-o: 

4.916454 

K-nc 

1.  S7Bf.fi  5 

?:-o. 

1. 315299 

k-02 

1.  100000 

K-OO 

2.  590449 

E-01 

2.823262 

E-Ul 

1. b4B9S5 

k-01 

S.27ED26 

Is. -03 

1.200000 

k-00 

2.40G352 

E-01 

-5.957774 

E-00 

1.341215. 

.E-01 

1.  165493 

k-02 

1.30(1000 

K-00 

1.  889666 

k-Ol 

-3.386384 

E-00 

1,064485 

K-01 

1.  183559 

E-02 

1. -100000 

E-00 

1.445752 

E-01 

1.043906 

E-01 

8.252229 

K-02 

8.  159877 

E-03 

1.500C03 

K-00 

1,375517 

E-Cl 

-2.B08S20 

E-00 

G.  22  73  6 6 

E-02 

3.3BB179 

E-03 

..  *'30000 

L-00 

7.  73931  1 

E -02 

3.21007B 

E-00 

4.511455 

E-02 

3.  554770 

E-0J 

1.  700000 

E-00 

5.42030G 

E-0  2 

-7.900078 

E-00 

3. 134969 

E-02 

4.824774 

E-03 

1.800000 

E-00 

3. 748272 

E-0S 

5.  173332 

E-00 

2. 117924 

E-02 

4.383155 

E-03 

1.000000 

k-00 

2.502468 

E-02 

-8.205769 

E-01 

1.445394 

E-02 

2.311fi07 

E-03 

2.000000 

k-00 

1. 7B7807 

E-02 

2.747514 

E-00 

1. 010137 

E-02 

7.  933B70 

E-04 

2.  :C030C 

1*  -30 

1. 173791 

E-02 

3. o3B513 

E-00 

0.  7994  99 

E-03 

2. 135224 

K-D3 

2. 200300 

1 -00 

7.41D72ii 

L -03 

-B. 114787 

E-00 

4.269860 

E-03 

2.525229 

E-03 

2. (00000 

Y -00 

4.635326 

E-03 

-5.217757 

E-00 

2.564192 

E-03 

1.B85208 

E-03 

2. 400000 

I -00 

1.007434 

E-03 

4.870997 

E-00 

1.600111 

E-03 

5. 088279 

E-04 

2. 500000 

E-00 

1.  374182 

E-03 

8. 843470 

k-00 

1. 113705 

E-03 

B.  781242 

E-04 

2. 1-00000 

E-00 

1.  18708! 

E-03 

7.  73B858 

E-00 

7.34T5h8 

E-04 

1.261554 

E-C3 

2. 700C0D 

E-C0 

6. 1704GB 

E-04 

-1.548346 

E-01 

4.318422 

E-04 

7. 290920 

E-04 

2.800000 

E-00 

2.4B2553 

E-04 

-3. 693822 

E-01 

2.075203 

E-04 

4.316742 

E-04 

2. 900000 

E-00 

5.  465200 

E-05 

-7. 520300 

E-01 

b. 258467 

E-05 

1.014008 

E-03  ' 

(.000000 

k-00 

0.0 

-1.000000 

E-02 

0.0 

0.0 

N .n.lN>r  or  paints*  3 1 
Number  of  intervals  fi 

Number  of  points  per  interval*  5 5 5 5 5 6 


Table  3 and  in  Figs.  4 and  5 for  the  case  of  the  unscattered  data  in  the 
same  six-interval  configuration  indicates  that,  as  would  be  expected, 
the  results  of  the  radiance  curve  fit  and  the  emission  coefficient  calcu- 
lations were  more  accurate  for  the  case  of  the  unscattered  data.  This 
conclusion  is  supported  further  by  the  data  listed  in  Table  9,  which 
shows  the  percentage  error  between  the  analytic  emission  coefficient 
values  and  the  calculated  values  for  the  two  cases.  Of  the  31  values 
listed  for  each  case,  27  are  smaller  in  magnitude  for  the  case  of  the 
unscattered  data. 
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Number  of  Points:  31 

Number  of  Intervals:  6 

Number  of  Points  per  I nterval:  5 5 5 5 5 6 

Data  and  Uncertainty 


1.20  1.80 
Displacement,  arbitrary  units 


Figure  8.  Intensity  for  scattered  e*2  test  data. 


Figure  9.  Emission  coefficient  profile  for  scattered 
b*2  test  data. 


,i J 

3. 00 
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Table  9.  Percentage  Error  Between  Analytic  Emission 
Coefficient  Values  and  Calculated  Values, 
Scattered  and  Unscattered  Data 


Displacement 


Percent  Error 
Scattered  Data 


Percent  Error 
Unscattered  Data 


0.0 

-1.407971 

E-01 

1.820905 

E-01 

0.1 

-1.  122049 

E-01 

1.603837 

E-01 

0.2 

-3.866298 

E-00 

1. 150176 

E-01 

0.3 

-4.325816 

E-00 

5.560210 

E-02 

0.4 

-7.  971595 

E-00 

4.809920 

E-02 

0.5 

4. 658263 

E-00 

7.776954 

E-02 

0.6 

-2.228545 

E-01 

7.425089 

E-02 

0.  7 

-3. 974722 

E-00 

5.  666415 

E-02 

0.8 

-5.884235 

E-00 

1.056823 

E-01 

0.9 

-7.074422 

E-00 

4.068959 

E-01 

1.0 

-4.447586 

E-00 

8.441329 

E-01 

1.  1 

-1.  729951 

E-00 

7.716219 

E-01 

1.2 

S. 380223 

E-01 

2.460397 

E-01 

1.3 

2.611250 

E-00 

-2.279747 

E-01 

1.4 

4.256720 

E-00 

2.937477 

E-01 

1.5 

5.217480 

E-00 

1.841442 

E-00 

1.6 

4.011329 

E-00 

1.755357 

E-00 

1.7 

6.  659431 

E-00 

-7.764359 

E-02 

1.8 

-3.355924 

E-00 

-2.300878 

E-00 

1.9 

8.567907 

E-00 

-1.494702 

E-00 

2.0 

-4.209002 

E-01 

3.922260 

E-00 

2.  1 

7.821717 

E-00 

6.823872 

E-00 

2.2 

-2.263196 

E-00 

5.385360 

E-00 

2.3 

-7.372242 

E-00 

9.468709 

E-01 

2.4 

-6.644632 

E-00 

-9.200700 

E-02 

2.5 

9.  127228 

E-01 

3.080837 

E-00 

2.6 

2.038233 

E-01 

1.  610981 

E-01 

2.  7 

2.396705 

E-01 

1.028519 

E-01 

2.8 

8. 656553 

E-00 

-1.447511 

E-01 

2.9 

-3.607077 

E-01 

-5.917934 

E-01 

3.0 

-1.000000 

E-02 

-1.000000 

E-02 
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3.3  TEST  DATA  WITH  DISCONTINUITY 


It  is  the  nature  of  least-squares  curve  fitting  techniques  to 
reduce  fluctuations  in  the  data.  Indeed,  the  techniques  are  designed 
to  possess  this  characteristic,  making  them  a highly  useful  tool  for 
reducing  the  effects  of  experimental  scatter  in  the  input  data.  However, 
when  a phenomenon  such  as  a sudden  change  in  slope  is  a vital  aspect  of 
the  raw  data  and  therefore  an  important  feature  of  the  emission  coeffi- 
cient values  to  be  calculated,  care  must  be  taken  in  manipulating  the 
parameters  to  achieve  a data  set  configuration  to  attain  optimum  accuracy 
in  the  resultant  emission  coefficient  values.  Since  the  basis  of  the 
inversion  technique  is  a polynomial  curve  fit,  it  is  reasonable  to  expect 
that  a good  radiance  curve  fit  will  be  difficult  to  achieve  when  the  data 
possess  an  abrupt  change  in  slope.  To  test  this  hypothesis,  a set  of  data 
was  constructed  by  using  the  following  form  for  the  emission  coefficient 
function: 


0.5 

, 0 < r < 

15 

25 

1 

15  < r < 

20 

\ -2 


(35) 


By  evaluating  Eq.  (1)  for  the  above  functional  form  of  e(r),  one  finds  the 
associated  analytical  radiance  function  to  be 


I(x) 


2 

( 400  - x ) 


50  - 1 

— sec 
x 


0 < x < 15 


15  < x < 20 


(36) 


Twenty-one  data  points  were  generated  for  input  data  using  Eq.  (36). 
These  input  data  are  listed  in  Table  10  along  with  the  associated  values 
of  the  emission  coefficient  found  by  evaluating  Eq.  (35)  at  the  appro- 
priate displacement  values.  Note  that  the  sharp  drop  in  magnitude  of 
the  radiance  between  the  displacement  values  of  15  and  16  is  reflected  in 
the  emission  coefficient  values  at  the  same  points.  The  initial  radiance 
data  are  shown  graphically  in  Fig.  10,  whereas  Fig.  11  displays  the  emission 
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coefficient  profile.  The  dotted  line  in  Fig.  11  represents  the  analytic 
results,  Eq.  (35).  The  smoothing  process  created  a curve  which  gives 
results  quite  different  from  the  analytical  results,  as  would  be  expected 
from  Fig.  10. 


Table  10.  Input  Data  Defining  a Curve  with  a Sudden 
Change  in  Slope 


Displacement 

Radiance  (Data) 

Analytic 

Emission 

Coefficient 

0 

20.00000 

0. 50000000 

1 

19. 97498 

0. 50000000 

2 

19.89975 

0. 50000000 

3 

19.  77372 

0. 50000000 

4 

19. 59592 

0. 50000000 

5 

19.  36492 

0. 50000000 

6 

19.  07878 

0. 50000000 

7 

18. 73499 

0. 50000000 

8 

18. 33030 

0. 50000000 

9 

17.  86057 

0. 50000000 

10 

17.  32051 

0. 50000000 

11 

16. 70329 

0. 50000000 

12 

16.00000 

0. 50000000 

13 

15. 19868 

0.50000000 

14 

14.28286 

0. 50000000 

15 

13.22876 

0. 50000000 

16 

2.01094 

0.09765626 

17 

1.63180 

0.08650519 

18 

1. 25285 

0. 07716049 

19 

0.83569 

0.06925208 

20 

0.00000 

0.06250000 
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Figure  10.  Intensity  for  twenty-one  point  sudden 
change  in  slope  data. 


1.80 

1.60 

1.40 

1.20 


Number  of  Points:  21 

Number  of  Intervals:  4 

Number  of  Poi  nts  per  I nterval : 5 5 6 5 

J Propagated  Uncertainty 


Figure  11.  Emission  coefficient  profile  for  twenty-six 
point  sudden  change  in  slope  test  data. 
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In  an  attempt  to  better  reproduce  the  steeply  sloping  portion  of 
the  data  between  the  displacement  values  of  15  and  20,  five  more  points 
within  the  interval  15  £ x ^ 20  were  generated  from  Eq.  (36).  This 
set  of  data  is  listed  in  Table  11  along  with  the  associated  analytic 
values  of  the  emission  coefficients.  The  initial  radiance  data  are 
shown  graphically  in  Fig.  12,  and  Fig.  13  displays  the  emission  coefficient 
profile  obtained  by  using  the  26  data  points  listed  in  Table  11  as  well  as 
the  analytic  results,  Eq.  (35).  A comparison  of  Figs.  11  and  13  shows 
that  the  sudden  change  in  slope  in  the  emission  coefficient  curve  is 
better  reproduced  using  the  second  set  of  data.  Nevertheless,  the  overall 
results  of  the  inversion  are  quite  unacceptable  and  seem  no  better  than 
the  results  illustrated  in  Fig.  11. 


Table  11.  Twenty-Six  Input  Data  Points  Defining  a Curve 
with  a Sudden  Change  in  Slope 


Displacement 

Radiance  (Data) 

Analytic 

Emission 

Coefficient 

0.0 

20. 00000 

0. 50000000 

1.0 

19.  97498 

0.  50000000 

2.0 

19. 89975 

0.  50000000 

3,0 

19.  77372 

0,  50000000 

4.0 

19. 59592 

0.50000000 

s.o 

19.36492 

0.  50000000 

6.0 

19.07878 

0. 50000000 

7.0 

18. 73499 

0. 50000000 

B.O 

18. 33030 

0. 50000000 

9.0 

17.8G057 

0.50000000 

10.0 

17.  32051 

0.50000000 

11.0 

16. 70329 

0. 50000000 

12.0 

16.00000 

0. 50000000 

13.0 

15. 19868 

0. 50000000 

14.0 

14.28286 

0. 50000000 

15.0 

13.22876 

0. 50000000 

15.5 

2.20671 

0. 10405830 

16.0 

2.01094 

0. 00765626 

16.5 

1.81998 

0.09182736 

17.0 

1.63180 

0.08650519 

17.5 

1.44389 

0.08163265 

1B.0 

1.25285 

0.0771G049 

18.5 

1.05341 

0.07304602 

19. 0 

0.83569 

0.06925208 

IB.  5 

0. 57455 

0.06574622 

20.0 

0.00000 

0. 06250000 
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Figure  12.  Intensity  for  twenty-one  point  sudden 
change  in  slope  data. 


Number  of  Points:  26 

Number  of  Intervals:  4 

Number  of  Points  per  i nterval:  8 8 5 5 


Figure  13.  Emission  coefficient  profile  for  twenty-six 
point  sudden  change  in  slope  test  data. 
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It  should  be  noted  from  examination  of  Figs.  10  and  12  that  the  por- 
tion of  the  radiance  curve  possessing  the  steep  slope  Is  not  veil  repro- 
duced in  either  case.  In  fact,  portions  of  the  fitted  curve  in  the 
steep  slope  region  lie  entirely  outside  the  standard  deviation  error 
bounds  on  the  radiance  data. 

It  is  clear  from  these  examples  that  it  is  not  always  possible  to 
adequately  fit  the  data  by  the  use  of  the  least-squares  spline  fit 
technique  described  herein.  It  is  reasonable  to  expect  that  when  the 
data  are  not  well  fitted  the  resultant  emission  coefficients  will  not 
be  accurate.  This  conclusion  is  supported  by  examination  of  the  data 
in  Tables  10  and  11  and  in  Figs.  11  and  13,  which  deal  with  the  emission 
coefficients  for  the  radiance  data  of  Figs.  10  and  12.  Consequently, 
in  such  cases,  a different  approach  to  the  inversion  or  a different 
set  of  fitting  functions,  something  which  has  the  potential  of  modeling 
the  discontinuity,  should  be  used. 

3.4  APPLICATION  TO  EXPERIMENTAL  DATA 

As  an  illustration  of  the  application  of  the  inversion  technique 
to  typical  experimental  data,  inversion  results  for  the  radiance  profile 
of  a selected  spectral  line  from  an  argon  arcjet  are  presented.  The 
details  of  the  experiment  and  implications  of  the  results  are  included 
in  Ref.  16.  The  radiance  profile  is  shown  in  Fig.  14.  Data  from  both 
sides  of  the  centerline  are  included  to  indicate  the  symmetry  of  the 
plume  and  to  provide  additional  data  for  the  least-squares  curve  fit. 

The  bars  shown  in  Fig.  14  represent  a typical  two-standard-deviation  bound 
for  the  data,  and  the  curve  represents  the  results  of  the  least-squares 
spline  fit.  The  results  of  the  inversion  and  error  propagation  are 
shown  in  Fig.  15.  The  error  bars  represent  the  two-standard-deviation 
uncertainty  at  each  of  the  radii.  The  largest  standard  deviation,  that 
on  the  centerline,  represents  approximately  a 10-percent  uncertainty  in 
the  corresponding  emission  coefficient.  The  results  represent  the  apparent 
best  fit  to  the  data,  as  determined  by  the  propagated  error,  and  represent 
physically  acceptable  results. 
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Distance  from  Plume  Centerline,  cm 


Figure  14.  Typical  lateral  radiance  scan  from  argon 
arcjet  at  X/D  = 2,  415.8  nm. 


Radius,  cm 


Figure  15.  Typical  emission  coefficient  profile  from  argon 
arcjet  at  X/D  = 2,  415.8  nm. 
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4.0  SUMMARY 

A method  of  performing  Abel  inversions  and  a method  of  determining 
the  propagation  of  the  associated  experimental  errors  has  been  presented. 
The  Abel  inversion  is  applied  to  the  problem  of  determining  the  radial 
distribution  of  emission  coefficients  from  observations  of  the  radiance 
from  a cylindrically  symmetric  radiating  source.  The  particular  scheme 
for  solving  the  Abel  inversion  problem  is  a least-squares  polynomial 
spline  fit  technique.  The  spline  fit  technique  involves  the  division  of 
the  raw  data  into  several  intervals  and  the  least-squares  curve  fitting 
of  the  data  points  in  each  interval  to  a sixth-degree  even  polynomial. 

The  ordinates,  slopes,  and  second  derivatives  of  the  polynomials  are 
required  to  be  continuous  at  the  interval  boundaries.  Thus  the  polynomials 
are  constrained  so  that  the  total  data  profile  is  smooth  and  yet  provides 
the  best  fit  of  all  the  data  in  the  least-squares  sense.  The  Abel 
inversion  of  the  resultant  polynomial  model  of  the  data  is  obtained 
analytically. 

, The  associated  error  propagation  analysis  is  developed  by  casting 
the  numerical  equations  selected  to  perform  the  curve  fit  and  integration 
into  a form  in  which  the  problem  can  be  viewed  as  a linear  transformation 
from  the  raw  data  to  the  inverted  results.  In  this  manner,  the  variance- 
covariance  matrix  of  the  raw  data  can  be  directly  transformed  to  the 
variance-covariance  matrix  of  the  emission  coefficients.  The  result  of 
the  error  propagation  analysis  provides  an  objective  basis  for  the  sub- 
jective determination  of  the  series  of  polynomials  providing  the  most 
nearly  correct  fit  to  the  raw  data  and  resultant  emission  coefficient. 

A computer  program  to  perform  the  least-squares  spline  fit  and  associated 
error  propagation  analysis  is  described  in  Appendix  A. 

To  determine  an  acceptable  least-squares  polynomial  spline  fit  for 
a particular  set  of  data,  there  are  generally  four  parameters  to  be 
considered:  (1)  the  total  number  of  data  points,  (2)  the  number  of 
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intervals  into  which  the  data  are  divided,  (3)  the  number  of  points  per 
interval,  and  (4)  the  displacement  distribution  of  the  data  points.  It 
has  been  shown  that  the  parameter  configurations  of  randomly  scattered 
data  which  yield  accurate  emission  coefficient  values  are  generally  the 
same  parameter  configurations  which  yield  accurate  emission  coefficients 
for  data  without  scatter.  However,  when  the  data  curve  possesses  an 
abrupt  change  in  slope,  it  is  generally  not  possible,  using  a polynomial 
function,  to  arrange  the  data  parameters  into  a configuration  that  yields 
an  acceptable  curve  fit. 

Although  the  development  of  the  error  propagation  method  has  been 
applied  specifically  to  a particular  least-squares  spline  fit  scheme,  the 
method  can  be  applied,  with  appropriate  modifications,  to  any  least- 
squares  polynomial  approximation  or  polynomial  spline  fit  technique. 
Therefore,  the  error  analysis  technique  can  serve  as  a means  of  comparing 
the  applicability  of  various  schemes  to  the  same  problem. 

It  should  be  noted  that  each  set  of  experimental  data  is  unique, 
and  the  choice  of  data  parameters  must  be  based  upon  an  analysis  of 
the  standard  deviations  generated  for  each  data  set.  The  error  analysis 
process  described  provides  an  objective  basis  upon  which  a subjective 
judgement  can  be  made  concerning  the  acceptability  of  the  results  obtained 
by  the  application  of  the  least-squares  polynomial  spline  fit  technique 
to  a particular  problem. 
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A.  1.0  GENERAL  INFORMATION 

The  analytic  and  numerical  approach  described  in  the  text  of  this 
report  has  been  coded  into  a computer  program  to  effect  the  Abel  inversion 
of  data  from  cylindrically  symmetric  sources  and  perform  the  associated 
propagation  of  experimental  errors.  The  purpose  of  this  appendix  is  to 
provide  the  description,  documentation,  and  user  manual  for  the  computer 
program.  The  program  described  herein  is  in  an  "as  developed"  state  and 
can  be  readily  modified,  as  required,  for  more  efficient  operation  for 
Abel  inversion  of  data  other  than  emission  data. 


A.1.1  DESCRIPTION  OF  PROBLEM 


The  physical  problem  is  the  determination  of  the  radial  distribution 
of  the  emission  coefficients  from  measurements  of  the  radiance  from  a 
cylindrically  symmetric,  optically  thin  radiating  source.  The  problem, 
illustrated  in  Fig.  A-1,  is  generally  expressed  mathematically  as 


)■(«).  2 /' 


f (r  ) r d r 

71  2~a 

l r — x ) 


where  y(x)  is  the  measured  radiance  as  a function  of  the  displacement 
x,  R is  the  overall  radius  of  the  source,  and  e(r)  is  the  radially 
dependent  emission  coefficient,  to  be  determined.  The  quantity  y(x) 
is  the  usual  experimental  measurement.  In  the  situation  described,  the 
emission  coefficient  can  be  expressed  as 


*(r) 


(dy/ dx  ) d 

71  TiT 

\ x — r ) 


x 
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However,  because  y(x)  is  subject  to  experimental  uncertainty,  there  is 
uncertainty  in  the  derivative  dy/dx.  Furthermore,  variations  in  dy/dx 
can  have  pronounced  effects  upon  e(r)  if  the  evaluation  proceeds  directly 
using  raw  data.  Consequently,  a smoothing  process  for  the  data  coupled 
with  a means  of  determining  the  effects  of  propagating  the  experimental 
uncertainty  through  the  smoothing  and  subsequent  inversion  is  required. 


Cylindrically 

Symmetric 


Figure  A-1.  Illustration  of  the  physical  problem. 
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The  smoothing  can  be  accomplished  readily  by  least-squares 
techniques,  and,  with  the  proper  choice  of  fitting  function,  the 
determination  of  the  emission  coefficient  can  become  analytic.  The 
details  of  the  mathematical  development  of  the  equations  for  smoothing 
and  subsequent  emission  coefficient  and  error  propagation  are  given  in 
the  body  of  the  report.  The  problem  may  be  summarized  by  noting  that 
a set  of  data  points  {(x,  y)} , where  x is  the  independent  variable  and  y 
is  the  dependent  variable,  are  to  be  curve  fit  to  a series  of  sixth- 
degree  even  polynomials.  Each  polynomial  is  to  be  valid  over  a specific 
range  or  interval  of  the  data,  and  adjacent  polynomials  are  to  be  smooth 
in  both  the  function  and  its  first  derivative;  the  second  derivative  is 
to  be  continuous  at  the  interval  boundaries.  The  series  of  polynomials 
is  to  be  further  constrained  so  that  they  provide  the  best  least-squares 
fit  to  the  data.  The  mathematical  problem  is  solvable  by  Lagrange's 
undetermined  multipliers.  Subsequent  to  the  least-squares  curve 
fitting,  the  last  interval  is  fitted  with  a polynomial  to  insure  zero 
slope  and  ordinate  at  the  outer  edge  of  the  data. 


The  dependent  variable,  y,  at  any  point,  x,  is  thus  expressible  by 
the  equation 


/ \ 2 4 6 

yU)=anx  + a|2x  4-  ai3x  + aj4  x 


(A—  1 ) 


where  the  subscript  i denotes  the  it  interval  such  that 

Zj_!  < x < Z.  (A-2) 

where  the  Z's  are  the  interval  boundary  points.  By  writing  Eq.  (A—  1 ) 
for  each  of  the  dependent  variables  (radiances)  measured,  one  can  evolve  a 
system  of  equations  linear  in  the  unknown  coefficients,  a^ , ai2,  a^, 
and  a^.  The  introduction  of  the  constraints  provides  additional  equations 
also  linear  in  the  unknown  coefficients  and  Lagrange's  multipliers.  The 
mathematical  curve-fitting  problem  is  thus  expressible  in  matrix  vector 
notation  as 

BA  = C (A- 3) 
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where  the  matrix  A is  the  (7n  - 3)  by  1 matrix  containing  the  coefficients 
of  the  polynomials  and  the  Lagrangean  multipliers  and  n is  the  number  of 
intervals  into  which  the  data  have  been  divided;  B is  a (7n  - 3)  by 
(7n  - 3)  symmetric  matrix  of  functions  of  the  independent  variable,  x, 
and  the  interval  division  points,  z;  and  C is  a (7n  - 3)  by  1 matrix  con- 
taining functions  of  the  independent  and  dependent  variables.  The 
matrices  B and  C are  described  in  greater  detail  below.  The  objective 
of  the  curve  fitting  technique  is  to  obtain  the  solution  of  Eq.  (A-3)  for 
the  column  vector  A. 


The  matrix  B may  be  partitioned  thus: 


(A-4) 


where  R is  a 4n  by  4n  block  diagonal  matrix  and  N is  a 4n  by  (3n  - 3)  block 
band  matrix,  each  of  which  may  be  further  partitioned  thus: 


R = 


P 

n 


(A-5) 


Q1 

- Q!  Qo 


•N  = 


- Q2  Q3 
_q3 


(A-6) 


54 


AEDC-TR-76-1 63 


are  each  4 by  4 symmetric  matrices  which  are 
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The  matrix  C in  Eq.  (A-3)  is  expressed  as  the  product  of  two 
matrices, 

C = GY  (A-9) 


where  G is  a (7n  - 3)  by  p matrix  of  certain  powers  of  x (the  p- 
independent  data  points)  and  y is  a p by  1 matrix  of  the  p-dependent 
data  points.  The  matrix  G may  be  partitioned  into  a block  diagonal  matrix 

'S, 


’1 


0 


(A- 10) 


where  each  is  a 4 by  m^  matrix 


S.  = 

j 


mj_  1+  l 


m . i + l 


‘j-l 


Xmj-1+2 


Xmj_  1 + 2 


m.  i +2 


mj_  1+  mj 

4 

6 

j-  1 +mj 


(A— 1 1 ) 


With  each  of  the  elements  of  B,  G,  and  Y thus  defined  for  a given  set 
of  data  (x,y),  solution  to  Eq.  (A-3)  is  immediate,  yielding  a vector 
containing  the  coefficients  of  the  polymonials  satisfying  the  least- 
squares  spline  fit  criteria.  The  solution  can  be  expressed  as 

A = WY  (A-12) 


where 


W = BT  1 G 


(A- 13) 
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With  the  dependent  data  thus  expressed  as  the  evaluation  of  a 
polynomial,  the  emission  coefficient, 


e (r) 


_1_  f R (dy/dx)  ax 
77  J 2 2 JS 

( x — r ) 


(A-14) 


is  analytic.  Since  the  range  of  integration  is  over  x,  and  the  dependent 
variable,  y,  is  expressed  by  different  polynomials  over  different  ranges 
of  x,  the  integral,  Eq.  (A-14),  is  expressed  as  the  sum  of  integrals, 
each  valid  over  a different  range;  i.e.. 


.w-i  r 


dx  + 


(x2-r2^ 


11 

£ J 


dy/dx 

~2  O’ 

Z (x  - r ) 
i-1 


Or,  substituting  the  polynomials  for  the  dependent  variables, 


t(r)  = - 


.2  2 % 
(x  - r ) 


- r > T „ 4 , 2 „ 2^ 

L 2 a2.j-  1 + T"  “3’j-l  U + * 


-i  i : i 11  ,2  2 ,!4  r- 


(A- 16) 


4 ,2-2.  6 ,„4  .22  _4. 

+ agi  ( x + 2 r ) 4.  — a4i  (3x  +4xr  +8r) 


which  is  linear  in  the  polynomial  coefficients.  Thus,  the  results  of 
evaluating  the  emission  coefficient  at  several  independent  data  points 
can  be  expressed  in  matrix  vector  notation  as 


E = MA 


(A- 17) 


where  E is  the  column  vector  of  emission  coefficients  evaluated  at 
m values  of  r,  A is  the  coefficient  matrix  determined  as  the  solution 
to  the  least-squares  spline  problem,  and  M is  an  m by  (7n  - 3)  matrix 
defining  the  coefficients  of  the  respective  elements  of  A. 
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Substituting  for  A,  the  emission  coefficient  may  be  expressed  as 

K = \1  l.r 1 G 'i  CA— 18) 

which  expresses  the  emission  coefficient  as  the  result  of  a linear 
transformation  of  the  data,  y. 

With  the  emission  coefficient  described  as  the  result  of  a linear 
transformation  from  the  data  space,  it  is  an  easy  step  to  provide  the 
transformation  of  the  uncertainties  of  the  data  to  uncertainies  of  the 
emission  coefficient. 


Let 


F = MBT'g 


so  that 


E = FY 


Then 


[El 


= FI Y 1 F 

L J C V 


CA— 1 9) 
CA-20) 
(A-2 1 ) 


where  the  [ ] ^ symbol  is  used  to  describe  the  variance- covariance  matrix 
of  the  parameter  enclosed. 

The  objective  of  the  computer  code  presented  herein  is  to  calculate 

the  elements  of  the  respective  matrices  M,  B,  G,  and  Y from  input  data 

{<x,  y)}  so  that  the  coefficient  vector  A and  the  transformation  vector 

F are  determined.  Further,  for  input  [y]  , the  [El  values  are 

cv  cv 

determined. 

A.1.2  LIMITATIONS 

As  written,  the  computer  program  is  subject  to  the  following 
restrictions  and  limitations: 
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1.  There  can  be  no  more  than  51  data  points. 

2.  There  can  be  no  more  than  10  intervals. 

3.  There  must  be  at  least  one  point  per  interval. 

4.  There  must  be  at  least  two  intervals. 

5.  The  variance-covariance  matrix  of  the  raw  data  must 
be  diagonal. 

6.  The  input  data  must  be  read  in  by  increasing  displacement. 

A.2.0  PROGRAM  DESCRIPTION 

The  program  in  its  present  form  requires  approximately  150  K bytes 
of  core  on  the  IBM  370/165  computer,  is  composed  of  six  Fortran  subroutines 
or  functions  and  a main  program,  and  will  perform  the  computations  for  31 
data  points  distributed  into  four  intervals  in  about  1/2  sec.  The  time 
required  per  case,  of  course,  varies  according  to  the  number  of  points, 
intervals,  and  points  per  interval. 

A.2.1  SUBPROGRAM  DESCRIPTIONS 

A short  description  of  each  pertinent  routine  used  in  the  computer 
program  is  listed  below. 

MAIN  PROGRAM 


The  main  program  provides  the  logic  for  the  computer  program  to 
execute  multiple  cases,  the  proper  calling  logic  to  the  various  routines 
which  effect  the  data  input,  inversion,  and  errors  propagation  analysis, 
and  the  program  summary  output.  The  program  utilizes  two  output  units: 
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the  line  printer  and  logical  unit  8.  The  logical  unit  8 provides  an 
additional  temporary  storage  device  so  that  the  results  of  the  analysis 
may  be  used  in  subsequent  analysis  programs  by  job  stepping.  Multiple 
data  cases  are  accomplished  simply  by  putting  the  proper  logic  inside 
a Do-loop. 

SUBROUTINE  INPUT 


As  the  name  implies,  this  subroutine  provides  for  the  input  of 
all  data  on  logical  unit  5.  The  subroutine  includes  calibration 
calculations  to  provide  for  conversion  of  the  input  raw  data  units  to 
physical  units.  The  input  data  in  raw  and  calibrated  form  are  output 
on  logical  unit  6.  All  communication  with  this  subrountine  is  through 
COMMON. 

SUBROUTINE  INVERT 


This  subroutine  provides  the  logic  to  perform  the  least-squares 
spline  fitting  of  the  data  and  determination  of  the  coefficients  and 
transformation  matrix.  The  bulk  of  this  work  is  accomplished  in  SUB- 
ROUTINE COVCAL.  However,  subsequent  to  the  call  to  COVCAL,  the  subroutine 
recalculates  the  coefficients  for  the  last  interval  in  order  to  assure 
meeting  the  constraints  of  zero  slope  and  ordinate  at  the  outer  edge  of 
the  data.  The  final  set  of  calculated  coefficients  are  output  on  logical 
unit  6,  and  the  percentage  difference  between  the  input  data  and  the 
results  of  the  curve  fit  are  calculated.  All  communication  with  this 
subroutine  is  through  COMMON. 

SUBROUTINE  COVCAL 


This  subroutine  provides  for  the  determination  of  the  elements  of 
the  matrix  B [Eqs.  (A-3) , (A-4) , and  (A-5)  through  (A-8)J , the  matrix 
G [Eqs.  (A-9)  through  (A— 11)],  the  matrix  M [Eq.  (A— 1 7 ) J and  the  matrix 
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F [Eq.  (A— 19)].  With  these  matrices,  the  polynomial  coefficients,  the 
emission  coefficients,  and  the  propagated  variance-covariance  matrix 
are  immediate.  The  calculations  proceed  calculating  sequentially  as 
follows : 

1.  The  array  G [Eq.  (A-9)  through  (A— 11)],  identified  in  Fortran  as  XT 

2.  The  array  B and  B 1 [Eqs.  (A-3)  through  (A-8)] , identified  in 

Fortran  as  (XTX) 

3.  The  coefficients  array  A [(Eq.  (A— 12)],  identified  in 
Fortran  as  AV 

4.  The  array  M [(Eq.  (A— 1 7 ) J , identified  in  Fortran  as  XTX) 

5.  The  array  F [(Eq.  (A- 19)],  identified  in  Fortran  as  XT) 

6.  The  array  [e]cv  [(Eq . (A— 2 1 ) ] , identified  in  Fortran  as  VC. 

All  communication  with  this  subroutine  is  through  COMMON. 

SUBROUTINE  EMSCAL 


This  subroutine  provides  the  logic  for  calculation  of  the  emission 
coefficients  at  each  of  the  radial  positions,  Eq.  (A-16),  numerically 
equal  to  the  input  displacement  positions.  The  calculation  proceeds  by 
Do-loop,  and  tracking  of  the  point  with  respect  to  the  interval  boundary 
points  is  maintained  so  that  the  correct  coefficients  are  used  in  the 
integral  evaluation.  Communication  with  this  subroutine  is  through  the 
argument  list:  No.  of  points,  No.  of  intervals,  interval  endpoints, 

displacement  array,  curve  fit  csefficient  arrays,  and  resultant  emission 
coefficient  array. 
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FUNCTION  EMFUN 

This  double  precision  function  performs  the  numerical  evaluation 
of  the  Abel  integral  at  either  the  upper  or  lower  limit  when  the 
intensity  is  described  by  a four-term,  sixth-degree  even  polynomial, 

Eq.  (A-16).  Input  arguments  include  the  radius,  the  upper  weighted 
by  the  exponent  of  the  corresponding  independent  variable  (that  is, 

2a^,  6a^).  These  weighting  factors  arise  from  the  differentiation 

of  the  polynomial. 

SUBROUTINE  MATINV 

This  is  a standard  service  routine  which  obtains  the  inverse  of  a 
matrix  by  Crout  reduction  with  partial  pivoting. 

A.2.2  LIST  OF  FORTRAN  VARIABLES 

A Spontaneous  transition  probability  (sec  ^) 

A1(I)  Curve  fit  coefficient  for  the  Ith  interval,  Eq.  (A-3) 

A2(I)  Curve  fit  coefficient  for  the  1^  interval,  Eq.  (A-3) 

A3(I)  Curve  fit  coefficient  for  the  I**1  interval,  Eq.  (A-3) 

A4(I)  Curve  fit  coefficient  for  the  Ith  interval,  Eq.  (A-3) 

Cl  Recalculated  coefficient  for  last  interval: 

A1 (NTVL) 

C2  Recalculated  coefficient  for  last  interval: 

A2 (NTVL) 
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C3  Recalculated  coefficient  for  last  interval: 

A3(NTVL) 

C4  Recalculated  coefficient  for  last  interval: 

A4(NTVL) 


CAL2 


CAL  3 
CAL4 

CALS') 
CAL6  j 

D 

DEN 

DISP 

EMIS 

EN 

ENDPT 

ENTEN 


Calibration  constant  for  intensity  according  to: 

ENTEN  . = (ENTEN,  - ZERO)  *CAL2 

cal  input 

Calibration  constant  for  x according  to: 

DISP  ..  = (DISP,  „ - CAL3)  * CAL4 

cal  input 

Not  used 

Intermediate  computational  variable 

Intermediate  computational  variable 

Array  of  displacements  at  which  data  are  taken 

Vector  E of  emission  coefficients 

Energy  level  (cm  ^),  not  used 

Array  of  interval  endpoints 

Vector  Y of  radiances  at  the  x locations  in 
DISP 


F 


Intermediate  variable 


FI  Assignment  statement  subprograms  for  the  computation 

of  Eq.  (A- 16),  Subroutine  COVCAL 


63 


AEDC-TR 

F2 

F3 

FAC 

G 

HEAD 

I 

IBOT 

ITIHE 

ITOP 

J 

J1 

J2 

J24 

JC 

K 


76-163 


Assignment  statement  subprograms  for  the  computation 
of  Eq.  (A-16),  Subroutine  COVCAL 

Assignment  statement  subprograms  for  the  computation 
of  Eq.  (A-16),  Subroutine  COVCAL 

Fractional  value  of  input  radiance  data  which  is 
to  be  taken  as  the  standard  deviation  of  the  radiance 

Statistical  weight 

Alphanumeric  header  for  identification 

Index 

Index 

Code  to  indicate  whether  x locations  and  interval 
sizes  are  the  same  from  one  data  set  to  the  next 

NTVL  + 1 

Index 

Index 

Index 

Index 

Index 

Index 
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K1 

Index 

K2 

Index 

K3 

Index 

K4 

Index 

L 

Index 

LIN 

Line  Counter  index 

M2 

Index 

M3 

Index 

MD 

7*NTVL  - 3 = 7n  - 3 

MSV 

Intermediate  variable 

NPNTVL 

Number  of  points  in  the  I1'*1 

NPT 

Total  number  of  data  points 

NPT1 

Computational  variable 

NSETS 

Number  of  sets  of  data 

NTLO 

Computational  variable 

NTVL 

Number  of  intervals,  n 

NTVL4 

Computational  variable 

interval,  m^ 
P 
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NTVL41 

NTVL42 

PERROR 


PI 

R 


RO 

SD 

SUM 

T1 

T2 

T3 

VC 

WL 

WORK 


Computational  variable 
Computational  variable 

Percentage  error  between  curve  fit  radiances 
and  input  data  values 

3.14159265 

A particular  x value  passed  to  EMFUN  from 
subroutine  EMSCAL 

Overall  radius  of  source,  R 

o 

Array  of  intensity  data  standard  deviations 

Intermediate  computational  variable 

Computational  variable  used  to  compute  Eq.  (A-16) 

Computational  variable  used  to  compute  Eq.  (A-16) 

Computational  variable  used  to  compute  Eq.  (A-16) 

Emission  coefficient  covariance  matrix  [El 

1 Jcv 

Wavelength  (angstroms) , not  used 

Array  containing  the  squares  of  the  interval 
end  points 
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X Intermediate  computational  variable  in  subroutines 

INVERT  and  EMSCAL:  array  used  first  to  store  the 
elements  of  the  matrix  W,  and  next  the  elements 
of  the  matrix  F [y]  in  subroutine  COVCAL 

XT  Array  containing  the  elements  of  the  matrix  G,  and 

then  the  elements  of  the  matrix  F 

XTX  Array  containing  the  elements  of  the  matrix  B, 

next  the  elements  of  the  matrix  B~\  and  finally 
the  elements  of  the  matrix  M 

Y A particular  endpoint  passed  to  EMFUN  from  EMSCAL 

Y1  Intermediate  variable 

YIP  Array  of  curve  fit  intensity  values 

YSPLN  Array  containing  data  intensity  vaues 

Z1  Intermediate  variable 

Z2  Intermediate  variable 

Z12  Intermediate  variable 

Z13  Intermediate  variable 

Z21  Intermediate  variable 

Z22  Intermediate  variable 
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ZERO  Calibration  constant  for  radiance  according  to: 

ENTEN  * = (ENTEN.  - ZERO)  CAL 2 
cal  input 

ZSPLN  Array  containing  squared  x values 

A.3.0  INPUT/OUTPUT 

The  input  consists  of  certain  control  parameters  defining  the  logical 
arrangement  of  the  physical  data,  the  radiance  data  and  the  corresponding 
position  (measured  from  zero  on  th£  centerline)  and  standard  deviation 
estimate,  and  calibration  factors.  In  addition,  there  are  four  input 
parameters  which  are  not  used  in  the  inversion  but  are  passed  on  to  the 
output  data  unit  for  subsequent  use.  When  narrow  spectral  line  data  are 
inverted,  the  four  additional  parameters  can  be  used  for  the  atomic 
constants  characterizing  the  radiation.  The  output  consists  of  input 
data,  calibrated  data,  certain  intermediate  calculation  steps,  and  final 
results.  The  output  is  clearly  labeled. 

The  principal  input  physical  parameter  is  the  radiance  (radiated 
power  per  unit  area  per  unit  solid  angle).  The  principal  output  physical 
parameter  is  the  emission  coefficient  (radiated  power  per  unit  volume  per 
unit  solid  angle).  The  specific  units  of  the  calculated  emission  coeffi- 
cient will  be  consistent  with  the  units  chosen  for  the  input  radiance  and 
displacement.  There  is  no  internal  unit  conversion  provided  other  than 
with  the  calibration  factors. 

A.3.1  INPUT  DATA  CARDS 

Card  No. 

and  Format  Fortran  Variable  Description 

1 . (13)  NSETS  Number  of  sets  of 

data  to  be  radially 
inverted. 
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2.  (20A4)  HEAD  Header  card  to  provide 

means  of  identifying 
uniquely  each  set  of 
input  data. 

3.  (6E12.0)  WL,A,G,EN  Four  variables  not 

used  directly  in  the 
calculation  but  passed 
through  to  the  output 
unit  for  subsequent 
use.  For  narrow 
spectral  line 
emission  data  these 
may  be  wavelength, 
transition  proba- 
bility, statistical 
weight,  and  energy 
level  respectively. 

4.  (313)  NPT  Number  of  data  points. 

NTVL  Number  of  intervals. 

ITIME  = 0 if  only  the 

radiances  have 
changed  from  the 
previous  data  set. 
= 1 each  time  a new 
set  of  data  is 
run. 

5.  (2613)  (NPNTVL(I), 1=1, NTVL)  NTVL  values:  each 

value  is  the  number 
of  data  points  in 
the  corresponding 
interval. 
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6.  (6E12.0) 

■ 

DISP(I) 

The  displacement 

ENTEN(I) 

radiance  and  standard 

♦ 

n. 

SD  (I) 

deviation  data  from 

I = 1 , NPT 

the  centerline  to  the 
outer  edge,  two  sets 
per  card. 

n+1.  (6E12.0) 

ZERO 

Radiance  data  calibra- 

CAL2 

tion  according  to 
I = (ENTEN  - ZER0)*CAL2. 

CAL3 

Displacement  data 

CAL4 

calibration  according 
to  X = (DISP  - CAL 3) 
*CAL4. 

CAL5 

Not  used. 

CAL6 

Not  used. 

For  multiple  data  sets,  repeat  cards  2 through  n+1 . 

A.3.2  OUTPUT 

The  printed  output,  on  logical  unit  6,  consists  of  five  pages/case, 
and  the  identification  and  logic  of  the  output  are  self-evident.  The 
printed  output  consists  of  1)  input  data,  2)  calibration  data,  3)  curve 
fit  coefficients,  4)  emission  coefficient,  and  5)  propagated  errors. 

In  addition  to  the  printed  output,  the  various  quantities  are 
written  (formatted)  on  logical  unit  8 for  offline  storage  of  results  or 
use  by  subsequent  job  steps.  This  output  is  listed  as  follows: 

Record  No. 

and  Format  Output  List  Description 

1 . (13)  J Case  No. 
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(20A4) 

HEAD 

Case  identification. 

(4D20.13) 

WL,A,G,EN 

Four  variables 
passed  from  input 
to  this  output 
record;  not  used 
in  the  radial 
inversion. 

(2613) 

NPT,  NTVL, 

(NPNTVL , I = 1,  NTVL) 

No.  of  points. 

No.  of  intervals 
and  No.  of  points 
in  each  interval. 

(4D20.13) 

(A1(I),A2(I),A3(I), 
A4(I),I  = 1,  NTVL) 

The  curve  fit 
coefficients  for 
each  of  the  intervals 

(4D20.13) 

(ENDPT(I) ,1  - 1.NTVL+1) 

The  endpoints  of 
each  interval. 

(4D20.13) 

(DISP(I) , ENTEN(I) , I = 1 ,NPT) 

Ordered  pairs  of 
displacement  and 
radiance,  two  pairs 
to  each  record. 

(4D20.13) 

(PEER0R(I),I  = 1 ,NPT) 

Percentage  error 
between  input  and 
fitted  radiances. 

(4D20.13) 

(EMIS(I) , I = 1 ,NPT) 

Calculated  emission 
coefficients. 
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(SD(I) , I = 1,NPT) 


Standard  deviation  of 
input  radiances. 


10.  (4D20.13) 


coefficient 

deviations. 


11.  (4D20.13) 


(VCLI) ,1  = 1 ,NPT) 


Emission 

standard 
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II 

III 

IB 

3C 

II 

III 

III 

III 

1 

EE 

1 

II 

IU 

II 

in 

:b: 

131 

331 

II 

a; 

131 

III 

id: 

131 

III 

IBC 

I 

in 

HI 

III 

III 

III 

ii: 

:i: 

:c: 

HI 

III 

IK 

13 

3 

11 

:d 

3 

II 

III 

II 

3QII 

III 

III 

III 

IEE 

II 

limn 

101 

IE! 

m 

in 

3QI 

III 

!D3£I 

III 

19! 

231 

III 

III 

III 

III 

III 

III 

llll 

III 

III 

III 

II 

III 

III 

III 

III 

III 

III 

mu 

u 

it 

IEE 

II 

imm 

101 

IEI 

m 

m 

III 

IEI 

!DHI 

III 

Ill 

in 

HE 

HI 

Ul 

III 

III 

III 

llll 

III 

III 

III 

II 

III 

III 

III 

III 

III 

!IJ 

mu 

u 

ii 
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A. 4.0  FLOWCHARTS 


MAIN  PROGRAM 


73 


SUBROUTINE  INVERT 


5 


AEDC-TR-76-163 


AEDC-TR-76-163 


Start 


Store  in  XT 
the  elements 
of  the  (7n  - 3) 
by  p matrix  G. 

Compute 

W = B-1G. 
Store  in  X. 

Define  bij  , i = 1,  4n 

j = 1,  4n,  where 

B = [ bii  ] , a (7n  - 3)  by 
{7n  - 3)  matrix.  Store 
these  elements  in  XTX. 


Compute 
A = WY. 
Store  in  AV. 


JJ  = 1 


Define  the  remaining 
elements  of  B.  Store 
in  the  corresponding 
positions  of  XTX. 


D054 
I = 1, 
< n 


Call  subroutine 
MATINV  to 

compute  the  elements 
of  B-*-  which  are  then 
stored  in  XTX. 


A1(I)  « AV ( JJ) 

A2  (I)  = AV  (JJ  + 1) 
A3(I)  = AV ( JJ  + 2) 
A4(I)  = AV ( JJ  + 3) 


SUBROUTINE  COVCAL 


AEDCTR-76-163 

A \ 


SUBROUTINE  COVCAL,  CONCLUDED 


AEDC-TR-76-163 


Note:  In  this  subroutine,  the  subscript  j (denoting 

interval  endpoint)  is  one  digit  larger  than 
the  subscript  i in  Eq.  (A-2) . 


SUBROUTINE  EMSCAL 


78 


Evaluate  the  term 


used  in  Eq.  (A-16) 


a3  (*2  + 2r2)  + TE  a4  (3x4  + 4x2r2  + 8r4 )] 


FUNCTION  EMFUN 
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A.5.0  OUTPUT 


CHECKOUT  OAT*  USINO  OEXPl-X*X)  FOR  INPUT  OAT* 
INPUT  OAT* 

NUMBER  OF  POINTS*  31 . 

number  of  intervals*  a 

NUMBER  OF  POINTS  PER  INTEHVALI 7 9 8 8 


INPUT  DATA  ARRAY 


0.4 

1 .0000000  00 

l.QOOOQOD-02 

1.0000000-01 

2.0000000-01 

9.9005000-01 

9.6079000-01 

9.900000D-02 

9.6100000-02 

3.0000000*01 

9.139300D-01 

9.1400000*02 

4. 0000000-01 

6.5214000*01 

8.5200000-02 

5. OOOOOOD-O 1 

7.7880000-01 

7.7900000-02 

6.000000D-01 

6.976800D-01 

6.9800000-02 

7. OOOOOOD-Ol 

6.1263000-01 

6.130000D-02 

a.ooooooD-oi 

5.2729000-01 

5.270000D-02 

9.0000000-01 

4,4860000-01 

4.490000D-02 

1.0000000  00 

3.6788000-01 

1,6800000-05 

1.1000000  00 

2.9820000-01 

2.9800000-02 

1.2000000  00 

2.3693000-01 

2.3700000-02 

I. 3000000  00 

1.8453000-01 

1.8400000-02 

1.4000000  00 

1,4086000-01 

1.41 0000D— 02 

1.5000000  00 

9.1054000  10 

1 . 050000D-02 

1,6000000  00 

Jj730500D-02 

7.7000000-03 

l.TOOOOOD  00 

5.55T6000-02 

5.600000D-03 

1.8000000  00 

3,9164000-02 

3.900000D-03 

1.9000000  00 

2,7052000-02 

2.7000000-03 

2.00.0000Q.JU] 

1.8316000-02 

1.8000000-03 

2.1000000  00 

1.2155000-02 

1.2000000-03 

2.2000000  00 

7.9071000-03 

7.910000D-05  _ 

2.3000000  00 

5.0418000-03 

5.0400000-05 

2.4000000  00 

3.1511000-03 

3,1500000-05 

2.5000000  00 

1.9305000-03 

1.930000D-0S 

2,6000000  00 

1*1592000-03 

1 « 160000D-05 

2.7000000  00 

6.8233000-04 

6.8200000-06 

2.8000000  00 

3.9367000-04 

3.9400000-06 

2.9000000  00 

2.2260000-04 

2.2300000-06 

3,0000000  00  _ 

1.2341000-04 

1.2300000-06 

BEGINNING .POINT 
0.0  6 

<2F  EACH  INTERVAL! 
.5000000-01  1.4500000 

00  2.2500000  00 
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CHECKOUT  DATA  USING  DEXP(-X»X>  FOR  INPUT  DATA 
INPUT  ( X , INTENS I TY.STO  DEV)  ARRAY 


0.0 

1.0000000  00 

1.0000000-02 

1.0000000-01 

9.9005000-01 

9.9000000-02 

2.0000000-01 

9.6079000-01 

9.6100000-02 

3.0000000-01 

9.1393000-01 

9.1400000-02 

A.  OOOOOOO-n 

8. 5214000-01 

8,5200000-02 

5.0000000-01 

7.7880000-01 

7.7900000-02 

6.0000000-01 

6.9768000-01 

6.9800000-02 

7.0000000-01 

6.1263000-01 

6.1300000-02 

8.0000000-01 

5.2729000-01 

5.2700000-02 

9.0000000-01 

4.486000D-01 

4,4900000-02 

1.0000000  00 

3.6788000-01 

3,6800000-02 

1.1000000  00 

2.982000D-01 

2.9800000-02 

1.2000000  00 

2. 369300D-01 

2.370000D-02 

1.3000000  00 

1.8453000-01 

1 • 840000D-02 

1 .4000000  00 

1.4086000-01 

1 .41 OOOOD-02 

1.5000000  00 

9, 1054000  10 

1.0500000-02 

1.6000000  00 

7.7305000-02 

7.700000D-03 

1.7000000  00 

5.5576000-02 

5.6000000-03 

l.eoooooo  oo 

3.9164000-02 

3. 9000000-03 

1.9000000  00 

2.7052000-02 

2.7000000-03 

2.0000000  00 

1.8316000-02 

1.8000000-03 

2.1000000  00 

1 • 215500D-02 

1.2000000-03 

2.2000000  00 

7.9071000-03 

7.9100000-05 

2.3000000  00 

5.0418000-03 

5.0400000-05 

2.4000000  00 

3.1511 000-03 

3.1500000-05 

2.5000000  00 

1.9305000-03 

1.9300000-05 

2.6000000  00 

1.1592000-03 

1.1600000-05 

2.7000000  00 

6.8233000-04 

6.8200000-06 

2.8000000  00 

3.9367000-04 

3.9400000-06 

2.9000000  00 

2.2260000-04 

2.2300000-06 

3.000000D  00 

1.2341000-04 

1.2300000-06 
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CHECKOUT  OAT*  USING  DEXP(-X»X)  FOR  INPUT  OATA 
CUBIC  COEFFICIENTS 


INTERVAL  START 

INTERVAL  END 

A 1 

A2 

*3 

0.0 

6.5000000-01 

-2.3*10600 

08 

1.0775*20 

10 

-5.1726060 

10 

5.3915260 

10 

^SSOftOCOrJIl 

1_.*S0090D  00 

*.229*810 

09 

-2.0918690 

10 

2.3289580 

10 

-5.2686800 

09 

1.450000D  00 

2.2500000  00 

-6. 1568760 

10 

7.2967050 

10 

-2.136*760 

10 

1.8108820 

09 

. 2.2500950  OQ. 

.3*0.000500  00 

l, *811*30 

11 

-6,7097350 

10 

9. *2*8080 

09 

-*,2201320 

OB 
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CHECKOUT  OAT*  USING  OEXP(-X*X>  FOR  INPUT  OATA 
INVERSION  RESULTS 


sta. 

DISPLACEMENT 

INTENSITY (OATA) 

INTENSITY ( CALC ) 

PERCENT  ERROR 

EMISSION  co cr 

1 . 

0,0 

XtOOOOOOO  00 

-2,3410600  08 

-2.3410650  14_ 

-2*4190980  09 

2 

1.0000000-01 

9.9005U0D-0 1 

-1.3147050  08 

-1.3279180  10 

-2.2175710  09 

3 

2.0000000-01 

9,6079000-01 

1.1759950  08 

1.2239870  10 

-1.6842780  09 

4 

3.000000D-01 

9.1393000-01 

3.5600450  08 

3.8953140  10 

-1.0226440  09 

s 

4.0000000-01 

8.5214U0O-O1 

3.8661010  08 

4.5369320  10 

-5.3473470  08 

6 

5.0000000-01 

7.7880000-01 

6.9294760  07 

8.897632D  09 

-5.560432D  08 

7 

6,0000000-01 

6,9768000-01 

-5.4318400  08 

-7,7855760  10 

-1*3103660  09 

8 

7.0000000-01 

6.1263000-01 

-1.0487060  09 

-1.7118100  11 

-2.6485610  09 

9 

8.0000000-01 

5.2729000-01 

-1.0002240  09 

-1.8969150  11 

-3,7914630  09 

10 

9. OOOOOOD-O 1 

4.4860000-01 

-2.3436090  08 

-5.2242720  10 

-4.5463200  09 

11 

1.0000000  00 

3.6788000-01 

1.3316860  09 

3.6198920  11 

-4.7670140  J9 

12 

1,1000000  00 

2.9820000-01 

3.6823460  09 

1.2348580  12 

-4.329970D  09 

13 

1.2000000  00 

2,3693000-01 

6.6676380  09 

2.8141810  J2_ 

-3.1552990  09 

1 A 

1.3000000  00 

1.8453000-01 

9.9633420  09 

5.3993070  12 

•1 .2400860  09 

IS 

1.4000000  00 

1.4086000-01 

1.3027370  10 

9.2484520  12 

. 1.2776750  09 

16 

1 .5000000  00 

9.1054000  10 

1.5075070  10- 

-8.3443810  01 

3.8959040  09 

IT 

1.6000000  00 

7,7305000-02 

1.5592340  10 

2,0169900  13 

5.8*02610  09 

18 

1.7000000  00 

5.5576000-02 

1.4575680  10 

2.622658D  13 

7.0*27750  09 

19 

1 .8000000  oo 

3.9164000-02 

1.2157890  10 

3*1043540  li 

7.4197050  09 

20 

1.9000000  00 

2.7052000-02 

8.6091330  09 

3.1824390  13 

6.929723D  09 

21 

2.0000000  00 

1,8316000-02 

4,3597210  09 

2,3802800  13 

5,5947120_09 

22 

2.1000000  00 

1 .2155000-02 

2.4260170  07 

1.9959010  11 

3.5772040  09 

23 

2.2000000  00 

7,9071000-03 

-3,5729390  09 

-4,5186460  13 

1.2127750  09 

24 

2.3000000  00 

5.0418000-03 

-5.5590270  09 

-1.1025880  14 

•5,5*43010  08 

25 

2.4000000  00 

3,  1511000-03 

-6,3219160  09 

•2*0062570  14 

~lt6?7^B10  09 

26 

2.5000000  00 

1.9305000-03 

-6.1181490  09 

-3.1692040  14 

-2.2331690  09 

27 

2.6000000  00 

1,1592000-03 

-5.1392170  09 

-4.4334170  14 

-2.376671D  09 

28 

2.7000000  00 

6.8233000-04 

-3.6489990  09 

-5. 347851D  14 

-2. 0R4670D  09 

29 

2.8000000  00 

3. 9367000-04 

-1.9915050  09 

-5.0588170  14 

-1,4*57070  09 

30 

2.9000000  00 

2.2260000-04 

-5.9892690  08 

-2.6905970  14 

-6.4?7049D  08 

__  31  _ 

3,0000000  00 

1,2341 00D-04 

0.0 

-1,0000000  02 _ _ 

^0,0 ; 
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CHECKOUT  OAT*  USINO  UEXPI-X»X)  FOR  INPUT  DATA 
ERROR  ANALYSIS  RESULTSI  STANDARD  DEVIATION 


DISPLACEMENT 

INTENSITY(UATA) 

STD  OEV ( INTENSI TY1 

EMIS  COEF 

STD  DEVIEhIS  COEF) 

0.0 

1.0000000  00 

1.0000000-02 

-2.4190980 

09 

1.0070710-01 

1.0000000-01 

9.9005000-01 

9.9000000-02 

-2.2175710 

09 

9.2976170-02 

2.0000000-01 

9.6079000-01 

9.6100000-02 

-1.6842780 

09 

7.1923150-02 

3,0000000-01 

9.1393000-01 

9.  1400000-02 

-1 , 0226440 

09 

4.437996D-02 

4.0000000-01 

8.5214000-01 

8. 5200000-02 

-5.3473470 

08 

2.6046750-02 

5.0000000-01 

7.7880000-01 

7,7900000-02 

-5.5604320 

08 

3.1480350-02 

6.0000000-01 

6.9768000-01 

6.9800000-02 

-1.3303660 

09 

3.4906640-02 

7,0000000-01 

6.1263000-01 

6.1300000-02 

-2.6485610 

09 

2.7101710-02 

6.0000000-01 

5.272900D-01 

5.2700000-02 

-3.7914630 

09 

1.8681150-02 

9.0000000-01 

4.4860000-01 

4.4900000-02 

-4.546320D 

09 

1.1866540-02 

l .ooooooo  oo 

3.6788000-01 

3.6800000-02 

-4.7670140 

09 

8.4821120-03 

1.1000000  00 

2.9820000-01 

2.9800000-02 

-4.3299700 

09 

8.9802740-03 

1.2000000  00 

2.3693000-01 

2.3700000-02 

-3. 1 552990 

09 

1 • 01 2472D-02 

1.3000000  00 

1.8453000-01 

1.8400000-02 

-1.2400860 

09 

9.7313310-03 

1.4000000  00 

1 .4086000-01 

1.4100000-02 

1.2776750 

09 

7.3695770-03 

1.5000000  00 

9.1054000  10 

1.0500000-02 

3.8959040 

09 

4,0595180-03 

1.6000000  00 

7.7305000-02 

7.7000000-03 

5.8802610 

09 

2.0515810-03 

1.7000000  00 

5.5576000-02  

5.6000000-03 

7. 082775a  09 

2.6411470-03 

1.8000000  00 

3.9164000-02 

3.9000000-03 

7.4397050 

09 

3.7069940-03 

1,9000000  00 

2,7052000-02 

2,7000000-03 

6.9297230 

09 

4.1107960-03 

2.0000000  00 

1.8316000-02 

1.8000000-03 

5.5947120 

09 

3.6971160-03 

2.1000000  00 

1.2155000-02 

1.2000000-03 

3,5772040 

09 

2.5221840-03 

2.2000000  00 

7.9071000-03 

7.9100000-05 

1.212775D 

09 

8.3402730-04 

_ 2,300  0 000_00_ 

5*0418000-03 

_ 5.0400000-05 

-5.SQ4301D  08 

_9. 4468390-04 

2.4000000  00 

3.1511000-03 

3.1500000-05 

-1.6274810 

09 

1.7725340-03 

2.5000000  00 

1.9305000-03 

1 *9300000-05 

-2.2331690 

09 

1.7634560-01 

2.6000000  00 

1.1592000-03 

1.1600000-05 

-2.3766710 

09 

9.9184920-04 

2.7000000  00 

6.8233000-04 

6.8200000-06 

-2.0946700 

09 

3. 735650D-04 

2.8000000  00 

3.9367000-04 

3.9400000-06 

-1.4657070 

09 

1.9629610-03 

2.9000000  00 
3.0000000  00 

2.2260000-04  __ 
1.2341000-04 

2.2300000-06 

1.2300000-06 

-6.4270490  08 
0.0 

1.0928190-01  _ 

0.0 
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A.6.0  FORTRAN  LISTING 


LEVEL  21.7  ( JAN  73  ) 


OS/ 360  FORTRAN  H 


COMPILER  OPTIONS  - NAME-  MA IN » 0PT-02  »L INECNT -56 .SIZE-OOOOK. 

SOURCt.EBCDIC.NOLlST.NOOECK.LOAD.MAP.NOEOIT.ID.aREF 
C THIS  PROGRAM  PERFORMS  RADIAL  INVERSIONS  BY  LEAST  SQUARES  SPLINE 
c fitting  the  raw  data  and  then  inverting  the  raw  data 


ISN 

ISN 

ISN 


0002 

0003 

0004 


ISN  0005 


ISN  __0006 


ISN 

ISN 


0007 

fl-QQ8_ 


IMPLICIT  REALMS (A-H.o-Z)  > 

INTEGER  HEAO<20) 

_0IMENSI0N  NPNTVLUO)  .OISP  1 51 1 .ENTEN  ( SI)  . ENQPT  < 1 2)  ,A1  ( 10) . A2(  101  • 
1 A3 ( ] 0 > * A4 ( 1 0 > .EMIS(51).YCALC<51 » ,PEWR0R(5H • SO (51 ) .VC (51 .51 ) 
_COMMON_DI  SP_iENTEN,ENOPT  * RO , A l . A2 , A 3 , A4 , YCALC  , PERROR . EM  I S . HE  AD . 

1 WL*A,G.fN.  SO. VC 

COMMON  NPT.NPNTVL.NTVL 


JSN 

0009 

( 

ISN 

0010 

C 

C 

_ ISN 

0012 

C 

c 

ISN 

0013 

_c 

ISN 

0014 

LSN 

0015 

ISN 

0016 

ISN 

0018 

ISN 

0019 

ISN 

0020 

ISN 

0021 

28 

_ JL$N 

0022 

ISN 

0023 

ISN 

0024 

30 

ISN 

0025 

ISN 

0026 

ISN 

0027 

ISN 

0028 

ISN 

0029 

ISN 

0030 

ISN 

0031 

ISN 

0032 

ISN 

0033 

ISN 

0034 

ISN 

0035 

ISN 

0036 

ISN 

0037 

ISN 

0038 

ISN 

0039 

ISN 

0040 

ISN 

0041 

100 

ISN 

0042 

C 

ISN 

0043 

1000  t 

READ (5* 1 005) NSETS 
00  100  J-l.NStTS 


CALL_ INPUT 
IFCJ.E0.1>ITIME»1 


CALL  INVERT 


CALL  EMSCAL (NMT,NTVL.ENDPT«OISP  « A 1 <1 ) . A2 < 1 ) ♦ A3 ( 1) • A4 < 1 ) ,EMIS) 
LIN-50 

DO  30  I-l.NPT 


IF (LIN.LT.SO)  GO  TO  28 
WRITE (6. 1000) HE AD 
WRITE (6* 1003) 

LIN-6 

CONTINUE 

WRITEI6.1004I  I.DISPtl)  ,ENTEN(I)  .YCALCdl  ,PERROR(n,EMIS(II 

LIN-LINM 

CONTINUE 

WRITE  <B(2000) J 

WRITE (8*2001) HE AO 

WRITE (8*2002) WL»A»G*EN 

*RITC(B,2000>NPT.NTVLf (NPNTVL(I) .l-l.NTVO 
ITOP-NTVLM 

WRITE (8.2002) (Al (I) *A2 (I ) . A3 ( I > , A4 ( I ) * 1-1 *NTVL) 

WRITE (8*2002)  (ENDPT ( I ) * I -1 . I TOP) 

WRITE (8.2002) (DISP ( I ) .ENTEN ( I ) • I-l *NPT) 

WRITE (8.2002) (PERROR (I).I-l.NPT) 

WRITEI8.2002) (EMIS(I).I-I.NPT) 

WRITE (8.2002) ( SO ( I > .1*1 »NPT ) 

WRITE (8.2002) (VC(I.I) .I-l.NPT) 

WRITE(6.1000)HEAO 
WRI TE (6.3000) 

WRITE (6.3002) 

WRITE (6*3001) (OISP(I) .ENTEN ( I ) .SD ( I ) .EMI S ( I ) .VC(I.I) .I-l.NPT) 

CONTINUE 

RETURN 

FORMAT  (•1*.20A4) 


DATE  76,271/08.51.03 
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PAGE  002 


ISN  0044 

ISN  0045 

ISN  0046 

ISN  0047 _ 

ISN  0046 

ISN  0049 

ISN  0050 
ISN  0051 
ISN  005? 

ISN  0053 


1003  FORMAT (IX, ‘INVERSION  RESULTS* .///l X. 'STA. » .SX, 'DISPLACEMENT* , 

2 S**»INTENSITY(OATA) •.5X*,INTENSITy<CALC) ».5X, 'PERCENT  ERROR*  * 

_ 3 5X » 'EMISSION  COEF'I 

1004  FORM*!  ( 1 X ♦ I3»  1PSE1 9»6)  _ ... 

1005  FORMAT < 13) 

2000  FORMAT  <26131 

2001  FORMAT (20A4) 

2002  FORMAT (4020. 1JI 

3000  FORMAT  (1X» 'ERROR  ANALYSIS  RESULTS?  _S.TANOARQ  PEylALLQNJ ,///j 

3001  FORMAT (IK* IPSE  19. 6) 

3002  FORMAT (8K.«0ISPLACEMENT».6X, • INTENSI Tr (OATA) ' ,2X»  *STD  DEV ( I NTENSI T 
IV) (tSXt'EMIS  LOEF 1 *6X. »STO  OEV(EMlS  COEF)*) 

ENO  _ 
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LEVEL  21.7  ( JAN  >3  I OS/360  RwTRAN  H ' OATE  76.271/08.51.06 

COMPILER  OPTIONS  - NAME  * MAIN.OPT*02.LINECNT«58.SIZE«OOOOK. 

_ SOURCE. EBCOIC.NOLlST.NODECK.LOAD. MAP, NOEDIT.IO.XREF 

ISN  0002  SUBROUTINE  INPUT 

C 

e 

C THIS  SUBROUTINE  OOES  THE  INPUT  Of  DATA  FOR  ABEL  INVERSION 

c~ 

c 


ISN  0003 
ISN  0004 

IMPLICIT  REAL#8  < A-H.O-Z) 
iNTtGER  HEAD <201 

ISN  0005 

r 

OIMENSION  NPNTVLUO)  *OISP(511  .ENTEN(bl),ENDPT(12l  ,A1  (10)  »A2<|0) 
1 A3  CIO) ,A4(I0i .EMI SCSI) ,VCALC<51 J .PEWRORCSl) . SO <51 ) .VC  151 ,51 ) 

ISN  0006 

COMMON  0ISP.ENTEN.EN0PT.R0.Al.A2.A3.A4.rCALC.PERR0R.EMIS.HEAD. 
1 WL*A«G.fN*  SO. VC 

ISN  0007 

c 

COMMON  NPT.NPNTVL«NTVL 

ISN  0008 

c 

READ (5. 1 000) Ht AD 

ISN  0009 

PEAO<5.inO?> WL.A.G.EN 

ISN  0010 

HEAD (5. 1001 ) NPT.NTVL.IT I ME 

ISN  0011 

REAL)  (5.  100)  .ERR* 100)  (NPNTVL  II)  »I*1.NTVL) 

ISN  0012 

REAUC5.1002.ENR.101I (OISP ( I ) .ENTEN ( I) »SD ( I ) »I«1 .NPT) 

ISN  0013 

READ(5.1002>ZtRO.CAL2.CAL3.CAL4,CALS.CAL6 

ISN  0014 

If  (CAL2.FQ.0.U)CAL2*1.0 

ISN  0016 

If (CAL4.EO.O.O)CAL4«1,0 

ISN  0018 

wRI  IE (6.20001  HEAD 

ISN  0019 

WHITtlft.ZOOll 

ISN  0020 

WRITE (6.2002) NPT.NTVL 

ISN  0021 

WRITE (6.2003) CNPNTVLCI) .l*l.NTVL> 

ISN  0022 

WRITE  <6.2005) 

ISN  0023 

c 

WRITE (6.2006)  CO  I SP C I » . ENTEN (I) .SO  1 1 ) . I ■ 1 . NPT ) 

c 

c 

CONVERT  INPUT  SCALE  READINGS  TO  INTENSITIES 

ISN  0024 

00  1 I ■ 1 » NPT 

ISN  0025 

1 DISP(n«OA0S(UlSP<I)-CAL3)»CAL4 

ISN  0026 

ENDPT(1)«0.0 

JSN  0027 

J-l  . . _ . _ 

ISN  0028 

JS-NPNTVLU) 

ISN  0029 

00  10  I ■ 1 .NPT 

ISN  0030 

If(l.NE.JS)  GO  TO  5 

ISN  0032 

J*J*1 

ISN  0033 

If  CJ.NE.NTVL*!) JS«JS*NPNTVL< J) 

ISN  0035 

If (I.NE.NPT)GU  TO  4 

ISN  0037 

ENDPT(J)-DISPINPT) 

ISN  0038 

IF  < J.EO.NTVL* 1 ) 60  TO  5 

ISN  0040 

WRITE (6. 3000 ) HEAD 

ISN  0041 

RETURN 

ISN  0042 

4 CONTINUE 

ISN  0043 

ENOPT(J)»(D!SP<n«DlSP(I*l)>*0.5 

ISN  0044 

5 CONTINUE 

ISN  0045 

fAC-S0Cl)/ENTtN(I) 

ISN  0046 

CNTfcN  C I > « (ENTtN ( I ) -ZERO ) *CAL2 

ISN  0047 

SDCl)»f AC*ENTtN(I) 
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ISN  0048  10  CONTINUE 

ISN  0 049  *RI TE<6,2004)  (ENOPT<I»,l«l.NTVL) 

ISN  0050  RRITFI6.2000)mEAO 

ISN  0051  *«irC(6»2007l 

ISN  005?  *RI f E (6*2006 ) (D l SP ( I > »CNTEN < 1 1 »S0 ( T ) • 1 ■ 1 *NPT) 

ISN  0053  ENOPTCNTVL*l)«niSP(NPn 

ISN  005**  NO*UISP(NPT» 

ISN  0055  RETURN 

ISN  0056  100  *RITE(6*4000> 

ISN  0057  4000  FORhat  ( »1N1*« 'ERROR  IN  NPTNVL'I 

ISN  0050  101  «RIIC<0, 4001)1 

ISN  0059  4001  FORMAT (///1X«'1*'«I5) 

ISN  0060  RETURN 

C 

ISN  0061  1000  FORMAT (?0A4) 

C 

ISN  006?  1001  FORMAT (?6I 3) 

C 

ISN  0063  100?  FORMAT (6E1?. 01 

C 

ISN  0064  2000  FORMAT! M ' .?0A4I 

c 

ISN  0065  2001  FORMAT (IX* 'INPUT  DATA') 

C 

ISN  0066  2002  FORMAT (/1X*«NUM0ER  OF  POINTS*' *I3*/1X« ’NUMBER  OF  INTERVALS* •• I 3) 

c 

ISN  0067  2003  FORMAT ( 1 X • 'NUMBER  OF  POINTS  PER  I NTERV AL I • • 26 1 31 

C 

ISN  0060  2004  FORMAT ( / 1 X * • 0lG I NN I NG  POINT  OF  EACm  INTERVAL* ••/( 1 X * 1 P 1 OE 1 3*6) ) 

c 

ISN  0069  2005  FORMAT (//IX* • INPUT  DATA  ARRAY' 1 

C 

ISN  0070  2006  FORMAT (1X*3(1 PE  13,6*0X1 I 

C 

ISN  0071  2007  FORMATUX. 'INPUT  ( X*  INTENSI  TY*STO  OEVI  ARRAY',///) 

c 

ISN  0072  3000  FORMAT  (•  1 • ,?0A4, /lx, 'J  OOtS  NOT  MATCm  rITH  NTVLM 

c 

ISN  0073  ENO 
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L£VEL  ?1.7  < JaN  73  ) OS/366  “TOMTRAfllf  DATE  76.271/08.51.10 

COMP ILElT"B»»T IONS  - NAME  * MAIN.OPT«02.LINECNT»5«.SIZE»OOOOK. 

SOURCt  tEBCOI C (NOL 1ST ♦ NODECK .LOAD* MAP  *NOED IT  *ID*XREF 


ISN 

0002 

c 

SUBROUTINE  INVERT 

tSN  0003 

IMPLICIT  REAL'S IA-H,0-Z)  7 

ISN 

0004 

INTEGER  ME AO (201 

ISN 

0005 

c 

01  ME.  MSI  ON  NPN1VLI 10)  .OISPI51)  .ENTEN (61 > .ENOPT  1 1?)  , *1  ( 1 0 ) ( 1 0 1 • 

1 A3<10) ♦A4(10> .EMI  St  51) .VCALC151 J . PEKROR « 5 ]> . SO < 51 ) . VC  15 1 . 5 11 

ISN 

0006 

DIMENSION  YSPLNUO 1 10  ) «ZSRLN  110*10  ) «MORK(  111 

ISN 

0007 

COMMON  OiSP.ENTEN, ENOPT.PO.il  .A2,A3.A*.YCAl£.PEP#0R.TmIS.HEAD, 

ISN 

1 WL'AtGtEN.  SO. VC 

0008 

c 

ISN 

0009 

c 

CALL  COVCAL 

ISN 

0010 

Zl-ENOPT(NTVL)**? 

ISN 

0011 

Z2«EN0PT<NTVL*1)»»2 

ISN 

0012 

Z1Z*Z1*21 

ISN 

0013 

Z13-Z12*Z1 

ISN 

0014 

Z21-Z2*Z2 

ISN 

0015 

Z22-Z2 1*Z2 

ISN 

0016 

C1«A1 (NTVL-1) 

ISN 

0017 

C2*A2<NTVL-1) 

ISN 

0018 

C3-A3<NTVL-1) 

ISN 

0019 

C4*A4 (NTVL“1 ) 

ISN 

0020 

Y1»CI*Z1*(C2*Z1*(C3  ♦C4*Z1)) 

ISN 

0021 

YIP*  C2*  2.0*  C3*  71  0.0*  CA*Z12 

ISN 

0022 

DFN-Z1-Z2 

ISN 

0023 

C4»(Y1P-2.0*Y1/OEN) / <OEN»OEN) 

ISN 

0024 

C3*Y1/<OEN*OEN»-C4*(Z1*2.0*Z2) 

ISN 

0025 

C2«-Z2»(2.0*CJ*3.0*C4«Z2) 

ISN 

0026 

C1«-Z2»(C2*Z2#IC3*Z2«C4) ) 

ISN 

0027 

A 1 (NT VL ) *C 1 

ISN 

0028 

A2(NTVL)*C2 

ISN 

0029 

A3 (HTVL ) *C3 

ISN 

0030 

c 

AM  (NTVL»«C4 

OISPLAY  COEFFICIENTS  FOR  THE  INTENSITY  CURVES 

ISN 

0031 

*RITE(6.1000)MEAD 

ISN 

0032 

1000  FORMAT  < ' 1 ' »20A4) 

ISN 

003  J 

MR  1TE  (6.  1 00 1 1 

ISN 

0034 

1001  FORMAT un ••cubic  COEFFICIENTS’*///) 

ISN 

0035 

VRITE(6*1002) 

ISN 

0036 

100J  FORMAT ( 2X , * INIERYAL  START •, 2X .' lNTEPYAL  END • , 7X . < Al • , 1 3X . • A2' . I2X , 

l’A3' *llX,»A4»j 

ISN 

0037 

*RltE(6. 1005)  (ENOPTII)  .ENOPT  U»l)  .11111  .A2II)  .A3I1I  .A*(!)  , 

2 I *1 .NTVL 1 

ISN 

0038 

1005  FORMAT  < IP6E15»6) 

c 

CALCULATE  INTENSITIES  AND  PERCENTAGE  ERROR 

ISN 

0039 

J«1 

ISN 

0040 

00  40  I«I.NPT 

ISN 

0041 

A«DISP ( I ) **2 

ISN 

0042 

IF(A.GT.EN0PT1J*1»*«2) J«J*1 

ISN 

0044 

Cl-Al (J) 

ISN 

0045 

C2»A2(J) 
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ISN 

0046 

C3«A3(JI 

ISN 

0047 

C4aA4< J) 

ISN 

004* 

YC6LCI 1 1 «CI (C2  ♦ *• < C3*C4*X ) I 

ISN 

0049 

pERROR ( I 1 *0«0 

ISN 

OOSO 

IF(ENTENM)  .EW.O.OIGO  TO  40 

ISN 

0052 

TERROR (II  * 1 00 . 0» ( YC4LC ( 1 1 -ENTEN (1)1 /tNTEN ( I ) 

ISN 

0053 

40 

COnT INUE 

ISN 

0054 

41 

continue 

ISN 

OOS5 

WET URN 

ISN 

0056 

END 
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LEVEL  21.7  ( JAN  73  > 


OS/ 360  FORTRAN  H 


COMPILER  OPTIONS  - NAME*  MAIN,OPT*02.LINECNT«5B.SIZ£«OOOOK, 

SOURCt, EBCDIC. NOLIST, NODECK .LOAO. MAP, NOEDIT. 10, XREF 

ISN  0002  SUBROUTINE  ENSCAL (NPT.NTVL.ENDPT.0ISP.A1 .A2.A3.A4.EMIS) 

ISN  0003  IMPLICIT  REAL*8(A-M.O-Z> 

ISN  0006  DIMENSION  EN"  " ' * ‘ * 


ISN  0005 
ISN  0006 
ISN  0007 
ISN  0008 
ISN  0009 
ISN  0011 
ISN  0013 
ISN  0014 
ISN  0016 
ISN  0017 
. LSN  QOJ  B 


NDPT(1)«0ISP<1>  «A1  (1)  .A2in.A3(ll  »A4(1)  ,EMISU> 


I MENS 
J*1 

00  60  I«l.NPf 

EMism*o.o  

x-dispci) 

ini.EO.NPT)  l»0  TO  60 
Ir  CX.GT.ENOPT (J*l ) ) J«JM 

EMIS(I)«EMFUNI2.0*A2(J) .4. 0#A3 I J) .6.0*A4 { J> .X.ENDPT ( J*1 ) > 

IF ( J.EO.NTVL I 60  TO  60 

IBOT* J* 1 

00  50  K-IBOT.MTVL 

EMlSM)»EMIS<i)*EMrUNC2.0*A2(K)  ,4.0#A3<K»  *6.0*A4(K)_,X»ENDPT  (KM ) ) « 


2 EMFUN(2.0*A2(K) ,4.0*A3(K) »6.0*A4(K> .X.ENDPT(K) ) 

ISN  0019  _ 50  CONTINUE 

ISN  0020  60  CONTINUE 

ISKL0021  RETURN 

ISN  0022  END 


DATf  T6. 271/OB. 51. 13 
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LEVEL  21.7  C JAN  73  » 0S/360  FORTRAN  M DATE  76.271/08.51.16 

COMPILER  OPTIONS  - NAME*  MAIN.OPT*02.LINECNT«5B,SI7E«OOOOK. 

SOURCt, EBCDIC. NOLI  ST .NODECK.LOAD.MAP. NOEDIT. ID. XREF 
ISN  0002  DOUBLE  PRECISION  FUNCTION  EMFUN ( A2  * A J . A4.R . Y ) 

1SN  0003  IMPLICIT  REAL*8U-H.0-Z) 

ISN  0004  PI*3. 14159265 

ISN  0005  F«V*Y-R»R 

ISN  0006  IF(F.LE.l.OO-JO)  GO  TO  10 

ISN  0008  F*OSORT(FJ 

ISN  0009  T1*A2«F 

ISN  0010  T2«A3*<Y#r*2.0»R»R»*F/3.0 

ISN  0011  T 3»F • 1 3. 0* Y##A*A. 0*R*R*V*Y*8, 0*H*#4 1/15.0 

ISN  0012  T3*AA*T3 

ISN  0013  20  CONTINUE 

ISN  0014  T 1 » T 1 *T2*T3 

ISN  0015  tMFUN«-Tl/PI 

ISN  0016  RETURN 

ISN  0017  10  EMFUN«0.0 

ISN  0018  RETURN 

ISN  0019  END 
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LEVEL  21.7  < JAN  73  ) 


OS/ 360  FOwTRAN  H 


COMPILE#  OPTIONS  - NAME*  MAIN,OPT*02,LINECNT*5»»SIZE*OOOOK, 

..u  SOURCE .CBCOIC  tNOL  1ST .NODEC* .LOAD, MAP*  NOEDI T , 1 0 . XWFF 

ISN  0002  SUBROUTINE  COVCAL 


ISN  0003 
ISN  0004 
ISN  0005 


ISN  0006 

ISN  0007 
ISN  0008 

ISN  0009 
ISN  0010 
ISN  0011 


ISN  0012 
ISN  0013 
ISN  0014 
ISN  0015 
ISN  0016 


•THIS  SUBHOUTINE  PERFORMS  THE  ERROR  ANALYSIS  OPERATIONS, 


IMPLICIT  REAL*8(A-M,0-Z> 
iNTtGER  HEAD ( 20 ) 

DIMENSION  NPNfVL(10).l)|SP<51),ENTENfi»n.ENDPT(12)tAl  (10)  ,A?(|0)  , 

1 A3  < | 0)  • A4  ( 1 0 ; *£HIS  (_51 ) • YCALC  C51 ) ,PEMR0R<51)  *50  1511  .VC  (SI  , 51 1 , 

I XTX<67,67) .XT (67.51) ,AV(40) 

COMMON  0ISP.ENTEN,ENDPT.H0,A1.A2,A3,A4.YCALC.PERR0R.EMIS,HFA0, 

1 WLtA.G.EN.  SO. VC 

COMMON  NPT.NPWTVL.NTVL 
COMMON  /REX/  A (67,51) 

FI  (A.B) *(2./3. 14159265) -DSORT ( (A*A)-(#«B) ) 

F2(A,8)» (4./( J.«3. 14159265) >*(DSQRT ( lA»A)-(0»fl) ) )• ( (A* A) ♦ (2,*R*B) I 
F3(A.B)*(6./(15.#3. 14159265) ) • ( DSORT ( ( A* A )-(R*B) ))•( I3.*A*A*A*A)« 

1 (4.*A*A*8*B) ♦ (8.*B*8*B*0) ) 

NTVL4»4*NT VL 
NTVL41*NTVL4.1 
NT  VL42*NTVL4-2 
MO*7*NTVL-3 
NT  1 0* 1 0*NTVL 


ISN 

0016 

D03 J* 1 .NPT 

ISN 

0019 

3 

XT<I,J)«0.0 

ISN 

0020 

*1*1 

ISN 

0021 

*2*0 

ISN 

0022 

*4«*10 

ISN 

0023 

L*  1 

ISN 

0024 

D05J* 1 ,NPT 

ISN 

0025 

*■0 

ISN 

0026 

OON1.K1.K4 

ISN 

0027 

lriN.eo.oiooTujoo 

ISN 

0029 

iFioisPui.ea. o.on 

ISN 

0031 

XT1I,J)*01SP(J)**K 

ISN 

0032 

GOTU201 

ISN 

0033 

200 

XT  ( 1 . J) *1.0 

ISN 

0034 

GOTU20I 

ISN 

0035 

210 

XT ( 1 , J) *0.0 

ISN 

0036 

201 

*■*♦2 

ISN 

0037 

4 

CONT INUE 

ISN 

0038 

*2**2* 1 

ISN 

0039 

IFIKJ.NE.NPNTKLIIM 

ISN 

0041 

K 1 **4 ♦ 1 

_ISN 

0042 

K4** 1*3 

ISN 

0043 

*2*0 

UN  00** 

L*LM 

ISN 

0045 

5 CONTINUE 

c 

ISN 

0046 

006 1*1, MO 

ISN 

0047 

D06J*1 ,MD 

DATE  76.271/08.51.20 
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ISN 

0048 

6 

XTXI|,J>*0.0 

ISN 

0049 

N?*4 

ISN 

ooso 

Kl«l 

ISN 

0051 

K.l 

ISN 

0052 

K3>0 

ISN 

0053 

82*0 

ISN 

0054 

1*1 

ISN 

0055 

JC*I 

ISN 

0056 

*3*1 

ISN 

0057 

MSV*NPNTVLCK) 

ISN 

0058 

009J* 1 »NT 10 

ISN 

0059 

SUM*0.0 

ISN 

0060 

007L*K| »MSV 

ISN 

0061 

lFlK2.ro. 0IGOIO202 

ISN 

0063 

IF  ll)ISP(L).FO«0.0)GOTO7 

ISN 

0065 

SUN«SUN*0ISP<LI*»K2 

ISN 

0066 

G0TU7 

ISN 

0067 

202 

SUM*NPNTVL  < K ) 

ISN 

0068 

L*MSV 

ISN 

0069 

7 

CONT INUC 

ISN 

0070 

K3*«3*l 

ISN 

0071 

XTX(T.JC)*SUN 

ISN 

0072 

XTX ( JC* I ) *SUN 

ISN 

0073 

IFIK3.FO. H2)GOT0203 

ISN 

0075 

jc*jcm 

ISN 

0076 

K2*K2*2 

ISN 

0077 

G0T09 

ISN 

0078 

203 

1-1*1 

ISN 

0079 

JC*I 

ISN 

0080 

K3«0 

ISN 

0081 

«2*N2-1 

ISN 

0082 

lF(H2.EO.0»GOrO204 

ISN 

0084 

K2*K2-«H2-M3> 

ISN 

0085 

N3*N3*1 

ISN 

0086 

G0T09 

ISN 

0087 

204 

K 1 *NSV ♦ 1 

ISN 

0088 

K«K*  1 

ISN 

0089 

HSV«MSV*NPNTVL<K1 

ISN 

0090 

M2*4 

ISN 

0091 

K?bO 

ISN 

0092 

*3*1 

ISN 

0093 

9 

C. 

continue 

ISN 

0094 

J*1 

ISN 

0095 

K«2 

ISN 

0096 

001SI*NTVL41«MO«3 

ISN 

0097 

XTX  1 1 . J) *l./2. 

ISN 

0098 

KTXlJ«n*XTX(l.J) 

ISN 

0099 

XTK(|«J*4)*»XIX(1«J) 

ISN 

0100 

XTX IJ*4.I ) *X  T A ( I . J+4 ) 

ISN 

"101 

XTX(I.J*n*«l./2.)»£NOPTlK)*rNOPT  (K) 

ISN 

0102 

XTX(J«l.l)*XTX4IiJ*l) 

ISN 

0103 

XTXCI.J*5)*-XIX(I. JM) 

ISN 

0104 

XTXt  J*5*n*xTA<!.J*5) 

ISN 

0105 

XTX(ItJ*?)*(|./2.)*E NOPT IK  1 **4 

ISN 

0106 

XTX  U«?.I)*XTA  ( I «J*2I 
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ISN 

0107 

XTXII,J*6)*-Xrx{ltj*2> 

ISN 

0108 

XTXIJ*6.I>«XTX(I,J*6> 

ISN 

0109 

XTX 1 I • J*3) » ( 1 */2* ) •ENDPT  (X)  *«6 

ISN 

0110 

XTXUO*I)«XTX(I,JO> 

ISN 

0111 

XTX<I.J*7)»-XTXII.J*3> 

ISN 

on? 

XTX(J*7,I)«XTXCI,J*7J 

ISN 

0113 

xixi!M«J*n«fcNOPT<K! 

ISN 

on* 

XTXIJM#J*1>«ATX!I*1.J*1! 

ISN 

0115 

XTxCI*l.J*5>«-XTX«I*ltJ*D 

ISN 

0116 

XTX ( 1*1) aXTX  < I* 1 * J*5l 

ISN 

0117 

XTXI 1*1 *J*2) •ENDPT (K)»«3 

ISN 

0118 

XTX(  J*2tI*l)»XTXU*l  *J*2) 

ISN 

0119 

XTX  (1*1  * J*6»—XTX<  1*1,  J*2> 

ISN 

0120 

XTXU*6»I*1)«*TXU*1.J*6) 

ISN 

0121 

XTX ( 1*1. J*3)*J. •ENDPT (K)»«5 

ISN 

0122 

XTX(J*3.I*l>a*TXn*l.J*3) 

ISN 

0123 

XTX  (1*1*  J*T>—  XTX  ( 1*1*  J*3) 

ISN 

0128 

XtX(J.T,IM>«ATX(lM..I‘T> 

ISN 

0125 

XTX  1 1*2* J*1 ) ■! • 

ISN  0126 

XTX(J*1,1*2)*1, 

ISN 

0127 

XTX(I*2.J*5>  — 1. 

ISN 

0128 

XTX( J*5«l*2»aXTX(I*2.J*5) 

ISN 

0129 

XTX(I*2*J*2) ■•••ENDPT IK) •ENDPT (K) 

ISN 

0130 

XTXU»2.I»2I-«TX(W..)«2)  

ISN 

0131 

XTX  Cl  *2 « J*6)—XTX<I*2*J*2) 

I5N 

0132 

XTX1J*6»1*2)«XTX(I*2,J*6) 

ISN 

0133 

XTXCI*2.J*3)«15.»EN0PTCK)**4 

ISN 

0134 

XTX«J*3.I*2)»XTX(I*2.J*3) 

ISN 

0135 

XTXll*2.J*7)«-XTXCI*2tJ*3» 

ISN 

0136 

XTXC  J*7,I*2)«XTX(J*iiJ*7) 

ISN 

0137 

J«J*4 

ISN 

0138 

K*K  • 1 

ISN 

0139 

15  CONTINUE 

ISN 

0140 

CALL  MATINVIXTX.MD) 

ISN 

0141 

001 31 ■i*«b 

ISN 

0142 

D01 3j* 1 t NPT 

ISN 

0143 

13  XU'J)*0. 

ISN 

0144 

003QI»1 |M0 

ISN 

0145 

DO30J«l*NPT 

ISN 

0146 

DO30K*ltNO 

ISN 

0147 

30  X(I.JI*XTX(I.K>»XT(K.JI •XII.JI 

ISN 

0148 

D052I-1.NTVL4 

ISN 

0149 

52  AV ( I ) «0 • 0 

ISN 

0150 

D053I ■ 1 *NTVL4 

ISN 

01S1 

0053 J* 1 t NPT 

ISN 

0152 

53  AVmaX(l*J>*tNTEN(J)*AVf  1) 

ISN 

0153 

JJm  1 

ISN 

0154 

OOSMM.NTVL  _ _ 

ISN 

0155 

A1  (D-AVCJJ) 

ISN 

0156 

A2I1 )aAV(JJ«l) 

ISN 

0157 

A3(I)«AV<JJ*2) 

ISN 

0156 

A4(I»«AV(JJ*3) 
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PAGE  OOP 


ISN 

0159 

c 

54 

JJaJJ*4 

ISN 

0160 

0031I«1*NPT 

]SN 

0161  _ 

PP31 J«lfMO 

ISN 

0162 

r 

31 

XTX ( I • J) *0 • 0 

ISN 

0163 

4 

NPT1«NPT-1 

ISN 

0164 

JI«1 

ISN 

0166 

J2*2 

ISN 

0166 

0019I«1,NPT1 

ISN 

0167 

D-OISPCIl 

ISN 

0(66 

1F<D.GT.EN0PTIJ1MHG0T016 

ISN 

0170 

GOTO-17 

ISN 

0171 

16 

J1«J1M 

ISN 

0172 

J2*J2*4 

ISN 

0173 

17 

'TXU.J2*— ri(ENOPT<Jl«l)'OlSP(II> 

ISN 

0174 

XTXI1.J2*1>— *2(EN0PT<J1M>  tOISP(I)) 

ISN 

0175 

XTX(I.J2*2)«-7  3(EN0PT(J1*1I  tOISPdll 

ISN 

0176 

IEI J2.EQ.NTVL421G0T019 

ISN 

0176 

Kl-JlO 

ISN 

0179 

J24*J2*4 

ISN 

0180 

001«MC«J26*NTVL4f4 

ISN 

0181 

XTX ( I .8 ) *F 1 (ENOPTCKl) ♦OISPI Ill-Fl (ENOPTIKI *11 ,OISP(I) ) 

ISN 

0182 

XTX<I*KM>«F2<ENDPT(Kl)*0ISP<m-F2(EN0PT(Kla|»*0ISP(ln 

ISN 

0183 

ATX(I«K*2)aF3IEN0PT (81 1*01 SP 1 1 » 1 -F3 < ENOP T 1 K 1 ♦ 1 1 * 0 1 SP C I ) 1 

ISN 

0164 

18 

6 1 *6 1 ♦ 1 

ISN 

0185 

C 

19 

continue 

ISN 

0186 

00201*1 *NPT 

ISN 

0187 

0020 J* 1 ,NPT 

ISN 

0188 

20 

XT  < I . J) *0.0 

ISN 

0189 

00211*1. NPT 

ISN 

0190 

0021 J*1 .NPT 

ISN 

0191 

00218a) .HO 

ISN 

0192 

C 

21 

XT 1 1 . J I *XTX (I .8 ) *X 1 K » J ) • XT ( I . J) 

ISN 

0193 

00221*1 .NPT 

ISN 

0194 

0022 J* 1 .NPT 

ISN 

0195 

?2 

XTX  C J« I ) "XT ( I . J) 

ISN 

0196 

4 

00251*1 .NPT 

ISN 

0197 

0025J*1,NPT 

ISN 

0198 

25 

X<I.J)*0.0 

ISN 

0199 

00261*1 .NPT 

ISN 

0200 

0026J*1.NPT 

ISN 

0201 

Q 

26 

X(I»J|*XT(I.J) *SD(JI *50 (J) 

ISN 

0202 

00271*1 , NPT 

ISN 

0203 

0027J*l .NPT 

ISN 

0204 

27 

VC«I.J)»0.0 

ISN 

0205 

00281*1 .NPT 

ISN 

0206 

0028J*1 .NPT 

I5N 

0207 

00288*1 .NPT 

ISN 

0208 

28 

VCH.J|aX<I.8)»XTX  (8.JI  *VCU.J) 

ISN 

0209 

00291*1 , NPT 
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PAGE  005 

isn~o2io  29"vc<i*n»o$u«nvc<i.in 

I SN  0211 RETURN 

ISN  0212  ENO 
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LEVEL  21.7  I JAN  73  I 


05/360  FOMTRAN  H 


DATE  ?6,2?l/0ife.5l.2V 


ISN  0002 
ISN  0003 
ISN  0004 
ISN  0005 
ISN  0006 
ISN  0007 
ISN  0000 
ISN  0009 
ISN  0010 
ISN  0012 
ISN  0013 
ISN  0014 
ISN  0016 
ISN  0017 
ISN  0019 
ISN  0020 
ISN  0021 
ISN  0022 
ISN  0023 
ISN  0024 
ISN  002S 
ISN  0027 
ISN  0028 
ISN  0029 
ISN  0030 
ISN  0031 
ISN  0033 
ISN  0034 
ISN  003S 
ISN  0036 
ISN  0037 
ISN  0038 
ISN  0039 
ISN  0040 
ISN  0041 
ISN  0042 
ISN  0043 
ISN  0045 
ISN  0046 
ISN  0047 
ISN  0048 
ISN  0049 
ISN  O0S0 
ISN  0051 
ISN  0052 
ISN  0053 
ISN  0055 
ISN  0056 
ISN  0057 
ISN  0058 
ISN  0059 
ISN  (1060 
ISN  0061 
ISN  0062 


9 OPTIONS  - NAME*  MAIN*OPT*02»LINECNT*5b.SI2E*0000«* 

SOURCE .EBCDIC  *NOL I $T .NODECK  *LOAO , HAP. NOEOI T* ID. XREF 
SUBMOOT  I NE  MAHNV  (A.  NN1 

IMPLICIT  REAL*BIA-M.O-Z)  

0 1 MENS  I ON  A (6  7.67) .LOCATE (67.3) 

00  1 N* 1 ,NN 
1 LOCATE  < N. 3) ■ U 
00  14  N*1 .NN 
ANAA*  0.00*0 
DO  6 1*1 .NN 

IF  ( LOCATE (1*3) .FQ. 0 > GO  TO  2 
GO  TO  6 

7  00  b J* 1 . NN 

IF  C LOCATE  I J.3) .EO.O  ) GO  TO  3 
GO  TO  5 

3 IF (UARSCA « I , Jl) .GT.AMAX)  GO  TO  4 
GO  TO  5 

4 AMAAvOABS  (AliiJI) 

IROtfs  I 

JCOL*  J 

5 CONTINUE 

6 CONTINUE 

IF  ( AHAX.GT. 1.00-15  ) GO  TO  7 
GO  TO  18 

7 LOCATE  I N* 1 I ■ IROW 
LOCATE (N, 2)*  JCOL 
LOCATE  C JCOL*  31  * 1 

IF  < IR04.NE.JC0L  ) GO  TO  8 
GO  TO  10 

8 JO  9 J*  l . NN 
SWAM*  A ( TROW* J) 

A ( I MOW • J) • A ( JCOL  t Jl 

9 A ( JCOL , J) » SWAP 

in  pivot*  a(jcol»jcoi.) 

A ( JLOL • JCOL ) * 1.00*0 
DO  11  J*1.NN 

11  A(jCOL.J)*  A ( JCOL • J)  / PIVOT 
DO  14  |*1. NN 

IF  ( I.NE.JCOL  ) GO  TO  12 
30  TO  14 

12  F«  A(I.  )COL) 

All. JCOL)*  O.UO*0 
JO  13  J*1.NN 

13  A I J . J ) s A ( I * J ) - F *A  ( JCOL  * Jl 

14  CONTINUE 

JO  17  N* 1 . NN 
L*  NN-NM 

IF  ( LOCATE (L.l) .N£. LOCATE IL*2)  » GO  TO  15 
GO  TO  ]7 

15  I«0**  LOCATE (L.ll 
JCOL*  LOCATE(L.2) 

JO  16  K ■ 1 « NN 
Swam*  aik.IROw) 

A ( K » l ROW ) ■ A (6, JCOL) 

A l A. JCOL)*  Swap 

16  CONTINUE 
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ISN  0003  IT  CONTINUE 

ISN  0064 RETURN  

ISN  0065  10  PRINT  1000 

ISN  0066 RETURN 

ISN  0067  1000  FORMAT  <1H  , 10X.  15HS1NGULAR  MATRIX) 

15N  0060 END s 
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F88-LEVEL  LINKAGE  EDITOR  OPTIONS  SPECIFIED  LET«MAP 

DEFAULT  OPTION(S)  U5EO  - SIZE* ( 100352* 12280) 


HOOULE  MAP 


CONTROL  SECTION  ENTRY 


NAME 

ORIGIN 

LfcNGTH 

NAME 

LOCATION 

NAME 

LOCATION 

NAME 

LOCATION 

NAME 

LOCATION 

MAIN 

00 

702 

INPUT 

708 

..  790 

INVERT 

E98 

476 

EMSCAL 

1310 

410 

EMFUN 

1720 

27C 

COVC AL 

1 9 AO 

1094L 

MAT  I NV 

122F0 

83A 

IHCLSQRT* 

12R30 

ISM 

OSORT 

12B30 

XHCFOXPI* 

12C90 

140 

FOXPI# 

12C90 

IHCECOMH* 

12DE0 

F61 

IBCOM# 

12OE0 

FDI0C5A 

12E9C 

INTS*TCH_ 

13026 

IHCC0MM2* 

13048 

650 

SEOOASO 

140C0 

IHCFCVTH* 

143A8 

1105 

AOCONA 

143A0 

FCVAOUTP 

14452 

FCVLOUTP 

144E2 

FCVZOUTR 

1463A 

FCVIOUTP 

149EE 

FCVEOUTP 

14EF0 

FCVCOUTP 

1 5 1 0 A 

INT6SWCM 

1S3F3 

ihcefnth* 

15560 

542 

ARITH0 

15560 

AOJSMTCH 

158FC 

1 

— 

• - - — 

IHCEEIOS* 

15AA0 

F28 

FI0C5* 

15AA8 

FI0CS8EP 

ISAAC 

IHCFI0S2* 

169D0 

52k. 

ihcerrm  • 

16FOO 

50C 

ERRMON 

16F00 

IHCERRE 

16F18 

IHCUOPT  • 

174E0 

300 

IHCETRCH* 

177E0 

2BE 

IHCTRCH 

177E0 

ERRTRA 

177E8 

IHCUAT8L* 

17A70 

630 

S8LANKC0M 

1R0A8 

5020 

REX 

100C8 

6ACB 

. . 

ENTRY  ADDRESS  00 

TOTAL  LENGTH  24890 

• •••USERPROG  DOES  NOT  EXIST  RUT  HAS*  BEEN  AODED  TO  OATA  SET 
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NOMENCLATURE 


Matrix  of  curve  fit  coefficients 


ki 


Coefficient  of  ^ in  the  k6*1  interval 


Linear  transformation  matrix,  defined  in  Eqs.  (13) 
and  (14) 

Vector  representing  right-hand  side  of  least-squares 
equations,  defined  by  Eqs.  (11),  (13),  and  (19) 

Matrix  of  emission  coefficients 


[E] 


cv 


Variance-covariance  matrix  for  emission  coefficients 


Represents  matrix  product,  MW 


F1’Fk’Fn 


Functions  in  developing  least-squares  equations, 
defined  by  Eq.  (8) 


Matrix  of  powers  of  x,  Eqs.  (19),  (20),  (21) 

2 

Radiance,  generally  watts/cm  /sr,  but  for  illustrative 
purposes,  arbitrary  units 


i,j,k 


Dummy  subscripts 


M 


Matrix  of  coefficients  of  the  curve-fit  coefficients 
Eqs.  (23)  and  (24) 


m 


Number  of  points  at  which  emission  coefficient  is 
evaluated 
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™k 

N 


n 


P. 

J 


PkCx) 


P 


Qj 


R 


r 


S. 

3 


Sk 


W 


x 


Y 


[Y] 


cv 


fch 

Number  of  points  in  the  k interval 
Submatrix  of  B,  defined  in  Eqs.  (17)  and  (18) 

Number  of  intervals 

Submatrix  of  R,  defined  in  Eq.  (16) 

Polynomial  in  k1"*1  interval,  Eq.  (5) 

Number  of  points 

Submatrix  of  N,  defined  in  Eq.  (18) 

Submatrix  of  B,  defined  in  Eqs.  (15)  and  (16),  or 
outer  radius  of  emission  source 

Radial  position,  generally  cm,  but  for  illustrative 
purposes,  arbitrary  units 

Submatrix  of  G,  defined  in  Eq.  (21) 

Sum  of  squares  of  residuals,  defined  in  Eq.  (6) 

Matrix  product  B *G 

Displacement,  generally  cm,  but  for  illustrative 
purposes,  arbitrary  units 

Column  vector  of  radiance;  also  used  in  place  of  I in 
Appendix  A 

Variance-covariance  matrix  for  radiance 
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y 


z 


k 


e 


s 


'k,1 


pk,2 


Fk,  3 


Scalar  member  of  Y,  radiance 
Right  side  abcissa  of  k^  interval 

3 

Emission  coefficients,  generally  watts/cm  /sr,  but 
for  illustrative  purposes,  arbitrary  units 

Lagrange's  undermined  multiplier 

Vector  of  mean  values  of  radiances 

Vector  of  mean  values  of  emission  coefficients 

Expected  value  operator 

3.14159265 

Constraint  on  polynomial  at 

Constraint  on  first  derivative  at  Z. 

k 

Constraint  on  second  derivative  at  Z. 

k 
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